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CRITICAL MEASURES FOR VECTOR ENERGY: GLOBAL 
STRUCTURE OF TRAJECTORIES OF QUADRATIC 
DIFFERENTIALS 


ANDREI MARTfNEZ-FINKELSHTEIN AND GUILHERME L. F. SILVA 


Abstract. Saddle points of a vector logarithmic energy with a vector polyno¬ 
mial external field on the plane constitute the vector-valued critical measures , 
a notion that finds a natural motivation in several branches of analysis. We 
study in depth the case of measures jl = (mi,M 2?M3) when the mutual inter¬ 
action comprises both attracting and repelling forces. 

For arbitrary vector polynomial external fields we establish general struc¬ 
tural results about critical measures, such as their characterization in terms of 
an algebraic equation solved by an appropriate combination of their Cauchy 
transforms, and the symmetry properties (or the S'-properties) exhibited by 
such measures. In consequence, we conclude that vector-valued critical mea¬ 
sures are supported on a finite number of analytic arcs, that are trajectories of 
a quadratic differential globally defined on a three-sheeted Riemann surface. 
The complete description of the so-called critical graph for such a differential 
is the key to the construction of the critical measures. 

We illustrate these connections studying in depth a one-parameter family 
of critical measures under the action of a cubic external field. This choice 
is motivated by the asymptotic analysis of a family of (non-hermitian) mul¬ 
tiple orthogonal polynomials, that is subject of a forthcoming paper. Here 
we compute explicitly the Riemann surface and the corresponding quadratic 
differential, and analyze the dynamics of its critical graph as a function of the 
parameter, giving a detailed description of the occurring phase transitions. 
When projected back to the complex plane, this construction gives us the 
complete family of vector-valued critical measures, that in this context turn 
out to be vector-valued equilibrium measures. 


1. Introduction and statement of main results 

1.1. Historical background. Statistical systems of many particles have been ob¬ 
ject of intense analysis for a long time, both from the perspective of physics and 
mathematics. The study of a particular type of interacting particle systems, the 
so-called determinantal point processes (and their cousins, Pfaffian processes), has 
been especially fruitful in the past thirty years. This success can be explained both 
by the ubiquitous character and flexibility of these models (describing the eigenval¬ 
ues of several random matrix ensembles, non-intersecting diffusion paths, random 
growth models, random tilings, to mention a few) and the introduction of new tools, 
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intrinsically related with the analytic theory of orthogonal polynomials and their 
generalizations. 

A common feature of these models is that either the joint probability density, 
the correlation functions, the normalization constant or the generating function 
can be expressed as a determinant (resp., a Pfaffian), and the right selection of the 
functions appearing in these determinants unveils the integrable structure of the 
underlying processes. Another unified property of these models is the possibility to 
put them in the framework of the so-called log gases, where the particles behave as 
charges on one or two dimensional sets, subject to the logarithmic interaction. 


The best known case is the spectrum of some unitary matrix ensembles 28 29 


described in terms of classical families of orthogonal polynomials on the real line. 
Their zeros, all real and simple, satisfy an electrostatic model that goes all the 
way back to Stieltjes 76 , and solve a minimization problem for the associated 


(logarithmic) energy. In other words, the zeros provide an equilibrium configuration 
on the conducting real line. Similar situation occurs for some non-intersecting paths 
models or for the six-vertex model in statistical mechanics [24| , to mention some 
more examples. 

Another classical framework for orthogonal polynomials (orthogonality on the 
unit circle), developed in the seminal works of Szego 78 , is connected to the analysis 


of the Ising model 30 31 


Further immediate generalizations of the problems above oblige us to leave the 
real line and extend the notion of orthogonality and the associated log gas models 
to the complex plane. The so-called non-hermitian orthogonal polynomials appear 
naturally in approximation theory, as denominators of rational (e.g. Pade) approx- 
imants to analytic functions 40 63 , or equivalently, in the study of continued 


fractions. Electrostatic models on a conducting plane are also associated with the 


polynomial solutions of certain linear differential equations 58 60 


Recently, non-hermitian orthogonality found its way to areas with a much more 
“modern” flavor, playing the crucial role in the description of the rational solutions 
to Painleve equations 12 17 , in theoretical physics [l,[2 20 25 and in numerical 


analysis 27 


Observe that due to analyticity, there is a certain freedom in the choice of the in¬ 
tegration contour for the non-hermitian orthogonal polynomials, which means that 
the location of their zeros is a priori not clear. As it follows from the fundamental 
and Stahl 74 75] (see also the survey [4] ), the limiting zero 


works of Nuttall 64 


counting measure of these polynomials can be characterized in terms of certain max- 
min problems of the logarithmic potential theory, where the energy is minimized on 
each compact set, but these minima are maximized over an appropriately chosen 
class of contours. These questions (whose simplest example is the Chebotarev’s 
problem 55 65 about the selection of sets of minimal logarithmic capacity) are 
closely related to the free boundary problems in the theory of partial differential 
equations. 

Motivated by his study of convergence of Pade approximants to functions with 
branch points, Stahl 74 75 not only proved that the zeros of the associated non- 
hermitian orthogonal polynomials distribute asymptotically as the equilibrium mea¬ 
sure on the set of minimal capacity within the appropriate class (proving in passing 
their existence), but also characterized these sets in terms of a symmetry property, 
now known as the S-property. These ideas were considerably extended by Gonchar 
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and Rakhmanov [40| to include the case when the orthogonality weight depends on 
the degree of the polynomial. The electrostatic model features now a non-trivial 
external field (or background potential) that persists in the asymptotic regime; in 
consequence, the notion of the 5-property had to be extended to deal with this 
external field. 

The resulting Gonchar-Rakhmanov-Stahl (or GRS) theory allows to formulate 
conditional theorems of the form: “if there exists a compact set with the 5-property 
corresponding to the associated electrostatic model, then the weak limit of the 
zero counting measures of the related non-hermitian orthogonal polynomials is the 
equilibrium measure of this set”. However, the existence of such compact sets is 
a problem from the geometric function theory, and it usually has to be settled 
in each concrete situation, basically from scratch. There are very few genuinely 
general rigorous results in this sense, for instance the most recent one of Kuijlaars 
and Silva 53 settling down the existence for polynomial external fields on the 


plane, and also the contribution by Kamvissis and Rakhmanov 45 46 which deals 
with the Green potential instead of the logarithmic potential. 

The notion of the 5-property found also a natural interpretation in the light 
of the Deift-Zhou’s non-linear steepest descent method for the Riemann-Hilbert 
problems 28 . One of the key steps in the asymptotic analysis is the deformation 


of the contours, which transforms oscillatory behavior into a non-oscillatory plus 
exponentially decaying one. This deformation consists in “opening of lenses” along 
the level sets of certain functions (known as g-functions). It is precisely the 5- 
property of these sets which guarantees the exponential decay on all non-relevant 
contours. 

The available log gas models on the conducting plane suggest to extend the 
analysis of the equilibrium distributions (local energy minima on the prescribed 
sets) to the study of saddle points of the logarithmic energy on the plane. This 
yields the notion of critical measures , studied systematically in 
original ideas in the unpublished manuscript 


58 , following the 


66 . The advantage of dealing with 


saddle points of the energy is the possibility to avoid the complications of the free 
boundaries, as well as to use variational techniques that go back to the work of 
Lowner and Schiffer. 

(see also [53]) is the fact that the Cauchy 


One of the main conclusions of 58 


transform of any critical measure satisfies (a.e. on the plane) a quadratic equation, 
known as the spectral curve in the context of random matrix theory. An immediate 
consequence of this fact is that the support of such a measure is a finite union of 
analytic curves, all of them being trajectories of the same quadratic differential on 
the plane. 

Quadratic differentials are known to play an essential role in the study of confor¬ 


mal and quasi-conformal mappings, univalent functions, etcetera. Teichmiiller 79 


based on the fundamental work of Grotzsch, formulated a principle according to 
which quadratic differentials and their singularities are associated with the solutions 
of extremal problems in geometric function theory. Although he did not prove any 
general result supporting his principle, we should give credit to Teichmiiller for 
putting spotlight on quadratic differentials. Schiffer 71 72 introduced several vari¬ 


ational techniques in the geometric function theory, allowing to find differential 
equations for the boundaries of the extremal domains, which in turn give a more 
formal connection with quadratic differentials. See 43 56 67 for further details. 
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As it was mentioned above, the Riemann-Hilbert asymptotic analysis for 2x2 
matrix-valued functions characterizing the non-hermitian orthogonality provides 
another natural connection with quadratic differentials: the “right” level curves for 
the ^-functions, that usually have a quite complicated structure, are trajectories of 
a suitable quadratic differential. Thus, the problem of existence of the appropriate 
curves can be reduced to the question about the global structure of such trajectories. 
This is not a simplification (the description of the global structure can be a formi¬ 
dable task), but it opens the gate to other techniques from the geometric function 
theory. These ideas have been successfully applied by Bertola 15 , being closely 


related (in fact, in some sense equivalent) to the max-min approach mentioned be¬ 
fore, see also 18 . In the same spirit and in a much more explicit form, quadratic 


differentials also appeared in similar asymptotic problems in 9 H)j )l2]pi9| 50 ,[57 . 
We emphasize that in all these problems, the underlying quadratic differentials are 
defined on the complex plane. 

Returning to our motivation rooted in the interacting particle systems, recent 
works have showed that more general models, such as the hermitian matrix model 

np6p2 


with an external source 21 23 , the two-matrix model 


non-intersecting paths [47||48| or the normal matrix model 26 |52| require consider¬ 
ing log gases with more than one group of charges. In this case, standard orthogonal 
polynomials are not enough, and we must turn to the so-called multiple orthogonal 
polynomials and to the associated vector equilibrium problems 39 . The connec¬ 


tion of these two notions was put forward first by Gonchar and Rakhmanov 38|39 


motivated by the analytic theory of Hermite-Pade approximants. The reader in¬ 
terested in further details is advised to look into the monograph by Nikishin and 
Sorokin 63 , which features a nice introduction to vector equilibrium problems from 
the approximation theory’s perspective. 

Although the analytic theory of multiple orthogonal polynomial is still in its 
infancy, two somewhat extreme cases are fairly well understood (see e.g. [5|80 ). The 
first one, known as the Angelesco system, corresponds to orthogonality with respect 
to measures living on disjoint sets, whose asymptotics is described by a vector 
equilibrium for repelling systems of charges. The other one, called the Nikishin 
system, where on the contrary the orthogonality is defined on the same set for all 
measures, is governed by a system of mutually attracting charges. In both situations 
the size of the charged ensembles is fixed. From the recent works of Aptekarev [3j 


and Rakhmanov 68 it became apparent that the “intermediate” cases require more 


complex log-gases models for their asymptotic descriptions, mixing both attracting 
and repelling charges, as well as allowing for charges “flowing” from one ensemble 
to another, see Sections |1.2| and [2] below. 

In the context of the non-hermitian orthogonality, multiple orthogonal polynomi¬ 
als yield naturally the notion of the critical vector-valued, measures , defined again 
as saddle points (and thus, solutions of some variational problems) for the total 
logarithmic energy. 

In the same vein as the standard non-hermitian orthogonality can be character¬ 
ized by a 2 x 2 matrix valued Riemann-Hilbert problem 28|36 , the multiple orthog¬ 
onality can be associated with larger (A: x k, k > 3) matrix valued Riemann-Hilbert 
problems 37 . Some heuristics (that in some situations can be made rigorous using 
the WKB approximation arguments 16l]) indicates that the asymptotics of multiple 
orthogonal polynomials (and the critical vector-valued measures) can be described 
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using higher-order algebraic equations. This is indeed the case, as it follows from 
some previous works and from the results of this paper. 

Some instances of the critical vector-valued measures have appeared in literature, 
basically in the study of the limit zero distribution of some multiple orthogonal 
polynomials satisfying non-hermitian orthogonality conditions. The most studied 
cases are those presenting strong symmetry, so that the support of these measures 
can be easily derived. Other works, such as [6.7 , consider non-symmetric situations. 
Their basic approach is to start from a general form expected for a cubic equation 
that should govern the system, and then make a genus ansatz, hence reducing the 
free parameters of the equation to a minimum. Although a vector equilibrium is 
present, due to the explicit form of the deduced algebraic equation its detailed 
study is not needed and usually bypassed. Moreover, in such situations there is no 
external field involved, which is indeed a considerable simplification, even in the 
scalar case. 

One of our main motivations is to understand the structure of the critical vector¬ 
valued measures with harmonic external fields. We follow the approach initiated 
in 53)58 , and show that the variational method, applied to critical measures, allows 
to deduce the corresponding algebraic equation (spectral curve) and the associated 
trajectories of a quadratic differential, now on a compact Riemann surface instead 
of the plane, whose analysis will be the central part of this paper. 


1.2. Main results: general polynomial external fields. Assume we are given a 
vector of three non-negative measures ft = (fj, 1,^2, M3 ), compactly supported on the 
plane. For the interaction matrix A and vector of real-valued functions (“external 
fields”) $ = (< pi,(j) 2 ,(/> 3 ) 



where $ 1 , $ 2 , $3 are polynomial^] and denoting by 

!{»,”) ■= [f log 1 


dn(x)dv{y) 

JJ F - y\ 

the logarithmic energy of two compactly supported measures ^ and u, we consider 


the total energy functional of the form 39 

3 


iJ r. 

E(p) = = Y aj,kI(Hj,l*k) + Y / <t>3 d N- (1-2) 

j,k =1 j= 1 J 

Matrix A is symmetric, singular and positive-semidefinite. Its only non-zero 
eigenvalue (of order 2) is 3/2, and the corresponding eigenspace is given by the 
vectors (tq, V 2 ,Vz) T € C 1 * 3 (in what follows, (-) T stands for the transpose) satisfying 


iq — v 2 — v 3 = 0. 


1 Although our methodology is, generally speaking, valid also for general positive-semidefinite 

interaction matrices and when $>'. are rational, we restrict ourselves here to the given matrix A 

and the polynomial external fields; these choices avoid some technicalities related to cumbersome 
notation and the order or coincidence of the poles. 
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As it will be seen below, the external fields in (1.1) must satisfy the compatibility 
condition of <f>' = ($ 2 , $ 2 , $ 3 ) being an eigenvector for the eigenvalue 3/2, that 
is, 

^{z)-& 2 {z) = %{z), :ec. (1.3) 

We complete the description of the electrostatic model by introducing the con¬ 
straints on vector-valued measures under consideration: given a parameter a £ 
[ 0 , 1 ], we require that p = (p±, P2, P3) satisfies 

M + M = l, M + |/X 3 | = a, |/x 2 | - |/x 3 | = 1 - a. (1.4) 

Here \p\ denotes the total mass (total variation) of the measure p. 

It should be pointed out that a particular case of this model (with all = 0 and 
a = 1/2) appears in the work of Aptekarev [3| and Rakhmanov (68 in their study of 
Hermite-Pade approximants 6||8,70 . Notice that A in <o> contains submatrices 
that are the interaction matrices both for the Angelesco and for the Nikishin cases. 

Let us denote by mes 2 the plane Lebesgue measure on C. We are interested in 
the existence of critical vector-valued measures within a family A4 a of measures 
p = (hi, H 2 , P 3 ) satisfying (1.4) plus the following conditions: 


each fij, j = 1, 2, 3, is a non-negative Borel measure, supported on a com¬ 
pact set in C of zero plane Lebesgue measure, 

mes 2 (supp Hj) = 0 , j = 1 , 2 ,3, 


and such that the energy E(p), defined in (1.2), is finite, 
the set 

S a '■= (J (supp fij n supp p k ) 

l<j<k <3 

is finite, 

card S a < 00 . 


(1.5) 


( 1 . 6 ) 


The vector equilibrium problems deal with the minimizers of the energy functional 


(1.2) over a family of measures p living on prescribed (and fixed) subsets of C, see 
e.g. 13 [39 41 and the references therein. As it follows from [13], the equilibrium 
measure (global minimizer of E{p)), if it exists in the class Ai a , is unique. 

Here we consider a natural extension of the notion of vector equilibrium for the 
conducting plane: the critical vector-valued measures, defined as follows. For t £ C 
and h £ C 2 (C), denote by p* the pushforward measure of p induced by the variation 
of the plane z 1 —> h t (z) = z + th(z), z £ C. 

Definition 1.1. We say that for a £ [0, 1], p = (p±, p 2 , P 3 ) is an a-critical measure 
(or a saddle point of the energy E(-)) if p £ M a , and 


lim 

t —^ 0 


E(p t ) - E(p) 


= 0 , 


(1.7) 


for every function h £ C 2 ( C). 


Sometimes, when the value of a is irrelevant or is clear from the context, we will 
simplify the terminology speaking simply about critical measures. 

In order to formulate our first results we introduce some notation. Given a 
(scalar) non-negative Borel measure A, compactly supported on the plane, we de¬ 
note by 

log t- -dX(x) 

\x-z\ 
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its logarithmic potential, which is harmonic in C \ supp A and superharmonic in C. 
Additionally, the principal value of the Cauchy transform 

C\z) ■■= lim f d\(x) (1.8) 

e_> ’ 0 + J\z-x\>e x ~ z 

is analytic in C \ supp A, and 
8U X 

2 —— (z) = C x (z), zeC\suppA, (1.9) 

az 

with 

d_ _ 1 f d__ ,d_\ 
dz 2 \ dx % dy) 

The interaction matrix A in 0 satisfies the following identities: 

l 1 1 °\ ( 1 \ 

2A = B t B = 3/ 3 - bb T , B = -1 0 -1 , b = -1 , (1.10) 

\0 -1 1 J \-lJ 

where / 3 is the 3x3 identity matrix. For a vector of measures jl = {p-i, p, 2 , Us) G 
M a , define the vector of functions |* = (£ 1 ,£ 2 ,6a) T through the equality 

£ = jB Q$' + C^ (1 .11) 

where C = (C^ 1 , C^ 2 , C^ 3 ) T . Alternatively, the functions £ 1 ,^ 2,£3 are given ex¬ 
plicitly as 


600 

60 ) 



<T>' ( z ) 

-^-A+C^{ z ) + C^(z), 

O 

- _ C w( z ) -C^{z), 

O 


( 1 . 12 ) 


&(*) 


4. 

3 3 


C^(z) + C^(z). 


It was proved in [58] Section 5] (see also [53]) that the Cauchy transform of a 
scalar critical measure is a solution of a quadratic equation mes 2 -a.e. on C, which is a 
characterizing property of such measures. Our first result is a natural generalization 
of this property to the critical vector-valued measures: 


Theorem 1.2. Suppose that for a € [0,1], jl = (/^i, p. 2 ,/r 3 ) £ A4 a is an a-critical 
measure in the sense of Definition 0 Then there exist a polynomial R and a 
rational function D with poles of order at most 2 such that the functions £i,&, £3 
in (1.11) satisfy mes 2 ~a.e. on C the algebraic equation 

£ 3 - R(z)£ + D{z) = 0. 

The polynomial R is given by 


(1.13) 


r{z) = ! E a j,k&j(zWk( z ) - 


*i(s) - 


j,k =1 


3 = 1 ' 


x — z 


dp,j(x), z G C. (1.14) 


Moreover, each of the measures /ii,/i 2 ,/r 3 is supported on a finite union of ana¬ 
lytic arcs, and they are absolutely continuous with respect to the arclength measure 
of their supports. 
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Each pole p of D, if exists, belongs to the support of at least two of the measures 
Mi; M 2 ; M 3 ; an d f or each such a measure pk its density satisfies 


Jf[ 7 \( z ) = \ Z ~P\ "( c + o(l)), 


as z -A p along supp/ifc, 


(1.15) 


for some nonzero constant c (which may depend on p and k), where is the order 
of p as a pole of D. 

We also prove a converse to Theorem |1.2| 

Theorem 1.3. Suppose that for a £ [0,1], fi = (pi, p 2 , M 3 ) £ A4 a is such that 
the measures p±,p 2 ,P 3 are supported on a finite union of analytic arcs, and are 
absolutely continuous with respect to the arclength. 

If for some polynomials V\, V 2 , V 3 with 

V((z) + Vf(z) + Vf(z) =0, z £ C, (1.16) 

the functions 


Uz) = V({z) + C^(z) + C^(z), 
6 (z) = Vf(z)~C^(z)-C^(z), 
Uz) = V 3 l(z)-C^(z) + C^(z), 


z £ C \ (supp pi U supp P 2 ), 
z £ C \ (supp/ii U supp P 3 ), 
z £ C \ (supp p 2 U supp P 3 ). 


(1.17) 


satisfy an algebraic equation of the form (1.131 for a polynomial R and a rational 
function D, then p is a critical measure for the energy functional 0-0 and 
the external fields defined by 

$i(z) := Vi(z)-V 2 {z), $ 2 ( 2 ) := Vi(z) - V 3 (z), $ 3 ( 2 ;) := V 3 (z) - V 2 (z). (1.18) 


equalities (1.171 take the form ( 1 . 12 ). 


Notice that with the definition (1.18) and the compatibility condition (1.16) the 


Few remarks are in order. 

Theorems 11.21 and 11.31 show that there is a clear connection between critical 
measures and cubic equations of a specific form. The fact that the equation is 
cubic has to do with the rank of the interaction matrix A in 0 , which is 2 : from 
the proofs it is natural to expect that critical measures (with > 3 components) 
for an energy functional with an interaction matrix of rank r should be related to 
algebraic equations of degree r + 1 . 


Furthermore, the compatibility condition (1.16) has to do with the structure of 
the interaction matrix ( 1 . 1 ), and finds its expression also in the absence of the £ 2 


term in the algebraic equation (1.13). Equation (1.16) and the structure of (1.17) 
should be seen as normalization conditions, ensuring that the underlying energy 


functional is of the precise form (1.2). Different choices for these identities would 


lead to a different interaction matrix in ( 1 . 2 ). 


The notion of the critical vector-valued measures can be extended to accom¬ 
modate also log-rational external fields and domains with non-trivial boundaries 
(“hard edges”, in the terminology of random matrix theory). The variations of 
the energy in these situations must allow for fixed points, see [58j for the scalar 
case. Fixed points and singularities of the external fields will appear among pos¬ 
sible points of blow-up of the density of the critical measures, and also as possible 


poles of the coefficients R and D in (1.13). 


Finally, a reader familiar with the notion of the scalar equilibrium and scalar 
critical measures might be puzzled by the fact that we allow for a blow-up of the 
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densities of such measures in the case of the polynomial external fields. Indeed, in 
the scalar case all critical measures in such circumstances have uniformly bounded 
densities on the whole complex plane 53 . This is no longer the case in the vector 


case, as the following two examples illustrate: they show that the components of 
the critical vector-valued measures can create “artificial hard edges”, obstructing 
each other and presenting unbounded densities when it was not expected, i.e. not 
created by the external field or by the boundaries of the conducting domain. 


Example 1.4 (Repulsive or Angelesco interaction). Consider the algebraic equa¬ 
tion 

£ 3 - £ - - = o. 

Z 

It has simple branch points at ±3^3/2, and a double branch point at z = 0, and 
three solutions specified by their behavior as 2 —» oo: 

a ( ^=-i+o(i), &(*)= i+^+o Qr). ^)=- i +^+°(?)> 


with £ 1 ( 2 ) + £ 2 ( 2 ) + £ 3 (z) = 0. A direct calculation shows that there are two 
measures, Mi> living on [—3\/3/2, 0], and H 2 {x) = /ii(— x), such that |//i| = |/Z 2 1 = 
1/2, and 


C tAl (z) = l-&(z), z€<C\[-3\/3/2,0], 

C^(z) = -l-&(z), z € C\ [0,3V3/2]. 


Hence Theorem 


1.3 


shows that fl = (ji i,/r 2 ,0) G M 1/2 is a 1/2-critical vector¬ 
valued measure, according to Definition |l.l[ corresponding to the interaction matrix 
(1.1) and the polynomial external fields given by 

$i(z) =-2 =-$ 2 (z), $ 3 (z) = - 2 z. 


Notice however that the densities of both measures / ij, j = 1,2, blow up at z = 0 
as |z|“ 1/3 . 

It is worth mentioning that the same pair of measures ( 111 , 112 ) solves the Ange- 
lesco-type vector equilibrium problem of minimizing under the assumptions 

supp/ii C (-oo,0], supp/j 2 C [0, +00), |/Ui| = |/u 2 | = M3 = 0. 

This equilibrium problem does exhibit a hard edge at the origin, and the proof is 
based on the same type of energy variation, but leaving z = 0 as a fixed point. It 
should be pointed out that a similar problem (but on finite intervals) was analyzed 
first by Kalyagin in [44]. 


Example 1.5 (Attractive or Nikishin interaction). Let us consider now the fol¬ 
lowing cubic equation, which (up to a certain rescaling) have appeared already 


11 


= a (L19) 

As in the Example |l.4[ it has simple branch points at ±3\/3/2, and a double branch 
point at z = 0. 


Let (81 — 12Z 2 ) 1 / 2 denote the branch of the square root in C \ [— 3x/3/2, 3x/3/2] 
with positive boundary values on (— 3\/3/2, 3\/3/2) from the upper half plane, and 
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let 


T(*) := -1 + 


3(81 - 12z 2 ) 1/2 \ 
—2.z 2 ) 


1/3 


also have positive boundary values on the same interval from the upper half plane. 
Then T is holomorphic in C \ [—3-v/3/2, 3\/3/2], and 


£iO) = \' e ' 


i/3 r (z 


1 


O 


1 

~2 I as 2 ; 


e 7ri/3 Y(z) _ 

is holomorphic in C \ [0, 3\/3/2]. Moreover, Cardano’s formula shows that is a 
solution of ( 1.19| l; the other two solutions are given by 

& (z) = €i(~z) and 6 (z) = -£i(z) - &(z). 

For the selected branch of T, the boundary values of on (0, 3\/3/2) from 
the upper and lower half plane, respectively, satisfy 

w{x) := (£i,+ “ Ci,-) (%) >0, x £ (0, 3\/3/2). 

Z7 Tl 

Thus, d/i 2 (x) := w{x)dx is a positive unit Borel measure on [0, 3v3/2], and its 
Cauchy transform is 

C^{z)=b{z)- 1 -, ze C\[0,3>/3/2)]. 

Analogously, /13 (z) := ^ 2 (— z) is a positive unit Borel measure on [—3-\/3/2, 0], and 
C>»(z) = -Z 2 (z) + |, ze C\[-3>/3/2),0]. 


Theorem 1.3 shows that /2 = (0, // 2 /2, yu. 3 /2) £ M 1/2 is a 1/2-critical vector-valued 


measure, according to Definition m corresponding to the interaction matrix ( | 1 . 1 [ ) 
and the polynomial external fields given by 

$i(z) =0, ®z(z) = z = -® 3 (z). 

Notice that now the densities of both measures /ij , j = 2,3, blow up at z = 0 as 
|z|- 2 / 3 . 

The next example illustrates that in Theorem |1.3| we cannot discard the degen¬ 
erate cases when the coefficient D in (1.13) is identically zero, or other situations 
when (1.13) becomes reducible over the field of rational functions: 

Example 1.6. If /i is a scalar critical measure for the external field ReE, where 

( 1 . 20 ) 


V is a polynomial, then there exists a polynomial Q for which 53 58 

2 


If we set 


then 


(^C^(z) + = Q(z), mes 2 ~ a.e. z £ C. 

V 

Ml = M3 = 0, fJ,2 = A*, $2 = V, = -$ 3 = —, 

EQ?) =I^\^) + J R eVd/J, 


( 1 . 21 ) 


which shows that the vector-valued measure jl = ( 1 * 1 , 112 , 113 ) is a-critical for a = 0 . 
The function £ 2 i n (1-12) becomes zero, whereas £ 1 , £3 are analytic continuations of 
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each other across supp/i. From the proof of Theorem |1.2| in Section [2] we conclude 
that the corresponding algebraic equation (3.9) is then reducible and has the form 

( 1 . 22 ) 


cr - Q{z)) = o, 


where Q is the polynomial in (1.20). 


Reciprocally, if the Cauchy transform of a probability measure ^ satisfies an 
algebraic equation of the form (1.20) for some polynomials V and Q , then the 


construction just carried out combined with Theorem 1.3 gives us that /x is a scalar 
critical measure for the external field Re V. 

In very much the same spirit, if we now set 

1*2 = M 3 = AO Mi = 0, $i = 2$ 2 = 2d> 3 = V, (1-23) 

then the energy of /I* is also given by (1.21), showing again that /x* is a-critical, 


now for a = 1. The respective functions £ 1,^2 are analytic continuation of each 
other across supp /j, whereas £3 = 0 on C, and the corresponding algebraic equation 


is reducible to the form (1.22). 


The critical measure in (1.23) violates the assumption (1.6), but the conclusions 


of Theorem 1.2 still hold. However, condition (1.6) cannot be simply dropped, as 


the following example shows: 


Example 1.7. Let $ 3 , $ 2 be two arbitrary polynomials. For a = 1/2 and the 
polynomial vector external field (j> = ( 0 i, fa, 4>i ~ 4 > 2 ) T , <t>j — Re < J > 7 , satisfying 
(1.3), consider the vector-valued measure jl = (mi, M 2 , /J 3 ), with 


1 

Mi = M2 = 2 M, 


M 3 = 0. 


The total energy of jl is 


e (m) = 7 I(M, M)+o / ( 0 i + <h)dM 


Thus, if /i is a scalar probability critical measure in the external field 

2 

</> = g(</>l + fa)- 


then (1.7) for jl takes place. In other words, jl is 1/2-critical in the external field <j>. 
However, the statement of Theorem |1.2| does not hold: each of the corresponding 


functions £ 1 ,^ 2 , £3 in ( 1 . 12 ) satisfies an algebraic equation of degree 2 , but if (j) 1 7 ^ 
</> 2 , these equations cannot be combined into a single algebraic equation of degree 
3. To be more precise, calculations in the spirit of the proof of Theorem |1.2| in 
Section [2] below show that if Q is the polynomial in the right hand side of the 
equation d lo]) for n (with V' = ($i + $ 2 )/ 3 ), then 


£1 = Q, 


6 = -^1 + * i -® 2 ). 


C.3 = - 2 (?i +4>7)- 


Clearly, the only assumption not satisfied by the measure jj in this example is 


( 1 . 6 ). 


A consequence of Theorem |1.2| is that the set 

E a = (supp Mi U supp M 2 U supp /x 3 ) \ S Q 


(1.24) 


is a finite union of disjoint analytic arcs, each of them contained in exactly one of 
the supports supp Mi, supp M2, supp M3- 
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Any orientation of an arc E £ S Q defines its left and right sides, as well as the 
left (or +) and right (or —) continuous boundary values of each £j on E, that we 
denote by £j±. In what follows, we use extensively (and without further warning) 
the notation mod 3, so that £o = £ 3 , £4 = £1 and so forth. 


Since £j’s are the three branches of the cubic equation (1.13), assumption (1.6) 
implies that for any curve E C there exists an index j £ {1, 2,3} such that 

&+( z ) =&-(*)> £(j+i)±( z ) = £Cj-i)=f( z )> ( L25 ) 

Critical measures are intimately connected with vector equilibrium problems and 
the vector symmetry (S'-property), as it is evidenced by the next theorem. 

Theorem 1.8. Given j £ {1,2,3}, let E be an open analytic arc in supp fij \ S a , 
not containing any branch point of (1.13). There exists a constant l = l( E) £ M. 
for which both the Euler-Lagrange variational equation (equilibrium condition ) 


J2 a o,kU» k (z) 




= 1 , 


(1.26) 


k= 1 


and the S -property 
d 


dn A 


k =1 / \k—l j 


hold true for z £ E, where n± are the unit normal vectors to E, pointing in the 
opposite directions. 

If, for j = 1,2,3, the constants 1(E) =: lj above are independent of the con¬ 
nected component E of supp/ij, and additionally satisfy h — l 2 = h, then the 
critical vector-valued measure is the equilibrium measure for the vector of contours 
(supp /ii, supp H 2 , supp /. 13 ), see 13 Theorem 1.8]. Examples of critical but non¬ 
equilibrium measures in the scalar case are well-known (see e.g. |58;). Using the 


construction of Example [L 6 j in particular (1.23), we can easily build vector critical 
non-equilibrium measures. 

The S'-property is the natural generalization of the scalar S'-property, taking 
into account that in those arcs the potentials of the measures Pj-i and Pj+i are 
harmonic. 

Critical vector-valued measures ft are also closely connected with trajectories of 
a quadratic differential on the Riemann surface 1Z = 1Z\ U IZ 2 U TZa, of the algebraic 
equation (1.13), with 

TZi =C \ (supp//i U supp/x 2 ), 

TZ 2 = C \ (supp m U supp/z 3 ), (1-27) 

TZ 3 = C \ (supp /i 2 U supp/z 3 ). 


The sheets are connected across arcs in H, 
for E C 


according to the following rule: if 


the condition (1.25) holds, then we connect the sheets TZj+i and TZj-i 
through E crosswise, identifying the ±-side of E on TZj+i with its =F-side on TZj-\. 

As usual, we regard the solutions £ 1 , 62 , £3 in (1.12) as branches of the same 
meromorphic function £ : 7Z —> C, defined by (1.13), so that function £,• is the 


restriction of £ to the sheet 7 Zj of the Riemann surface 1Z. 
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Theorem 1.9. Let 

(£ 2 ( 2 )-&(z), on Tli, 

on iz 2 , 

[^ 1 ( 2 ) - 6 (-), onTZ z . 

Then Q 2 extends to a meromorphic function on TZ. Moreover, 

vo = —Q 2 (z)dz 2 


(1.28) 


is a meromorphic quadratic differential on TZ with possible poles only at the points 
at 00 . Finally, for j = 1,2,3, each arc o/supp/Xj is an arc of trajectory of vo. 

As it was mentioned in Section [l.l| the connection between quadratic differentials 
and critical (or equilibrium) measures is not new. For scalar critical measures, the 
quadratic differential can be globally projected on the complex plane, so it is not 
necessary to consider it on a non-trivial Riemann surface. Even for vector-valued 
measures, this connection has been exploited before, see for instance 51 54 , but 


only locally on the plane. To our knowledge, Theorem |1.9| is the first time this 
connection is stated on the whole surface 7 Z, which provides a powerful tool for 
construction of critical measures. 


Remark 1.10. Assume that at z = p we have (1.151 with v £ {1/3,2/3} (see 


Theorem 1.2 1 . Then p is a triple br anch point of 1Z , so that the local coordinate for 
7Z at z = p is 2 = p + w 3 . Theorem 1.9 asserts that supp Pj lives on trajectories of 


the quadratic differential vo defined in ( 1 . 281 , which in this local coordinate takes 


the form vo = w A 6 "(const +o(l))dw 2 , const 0 . 


If v = 1/3 (or equivalently if the polynomial D in (1.13) has a simple pole at p), 


then by the local structure of the trajectories, their canonical projection from 1Z 
onto C consists of two analytic arcs intersecting at p perpendicularly (a “cross”). 
As Example |1.4| shows, this configuration can correspond to the interaction of two 
repulsive measures, and p can potentially belong to the intersection of the three 
components pj. 

However, for v = 2/3 (which means that D has a double pole at p) the canon¬ 
ical projection of the trajectories from 7 Z onto C at p consists of a single analytic 


arc through this point. It implies, due to (1.6), that in this case p belongs to the 


intersection of at most two (and then, exactly two) components of /x; moreover, 
it will be the endpoint of an arc in each of the supports. A quick analysis of the 
cubic branching shows that these two positive measures cannot interact repulsively 
(“Angelesco” interaction), and hence, must attract each other (“Nikishin” interac¬ 
tion). Thus, Example |1.5| represents (at least, locally) the most generic case of the 
2/3-blowup of the density of a critical vector-valued measure. 


Remark 1.11. We should point out that (1.6) plays a substantial role for Theo¬ 
rem 1.2 but the conclusion of Theorem 1.3 still holds true even if we drop (1.6) 


completely. In this situation, if X is an analytic arc contained in the support of 
at least two of the measures p\, p 2 , P3, and p 2 7 ^ P3 on X, then we can choose an 
orientation on X for which 


£i+00 = ( 2 ), 6+(z) = £i-0), &+(z) = &-(*). 2 <EX. 


(1.29) 


The Euler-Lagrange equation (1.26) holds true for the three measures Pi,p 2 ,P3 


(and possibly different constants l), and with respect to the orientation (1.29) the 
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S'-property now becomes 





k= 1 


\k=l 



The conclusion of Theorem 1.9 is also valid in this situation (with appropriate 


gluing of the sheets of TZ across the overlaps of the supports). However, we do not 
know whether there are positive measures with non-trivial overlapping supports for 


which the functions ( 1 . 12 ) are the three solutions to the same cubic equation as in 
Theorem 11.31 


1.3. Main results: the cubic case. In order to illustrate the results stated above 
we analyze in depth the interesting and highly non-trivial cubic case, corresponding 


to the energy functional ( 1 . 1 ) ( 1 . 2 ), with 

®i(z) = $ 2 ( 2 ) = 2 3 , $ 3 (z) =0, z£C, 

so that the external fields are 

</>( z ) = <Ma> = 02(z) = Rez 3 , 0 3 (z)=O, z£ C; 


(1.30) 


(1.31) 


in what follows, we consider only the values of the parameter a £ [ 0 , 1 / 2 ); see the 
discussion in the Remark 11.151 and also in the introduction to Section |3] 

We build a continuous one-parameter family of critical vector-valued measures 


p by choosing in (1.13) the coefficients 


R{z) = 3 z 4 — 3 2 — c, D(z) = —2 z 6 + 3z 3 + cz 2 — 3a(l — a), 
where c = c(a) is the real parameter given by 

c = -(^ ( 1 - 4 a ( 1 -a )) 2X3 


(1.32) 


(1.33) 


Theorem 1.12. For a £ [0,1/2), there exists a one-parameter family p = ft a £ 
A4 a of critical measures for the potentials (1.31) for which the corresponding alge¬ 


braic equation (1.13) has coefficients (1.32). 


Moreover, for any such a critical measure, if the associated algebraic equation 


(1.13) defines a Riemann surface of genus 0, then its coefficients are given by (1.32) 


for some choice of the cubic root in (1.33). 


Remark 1.13. As it follows from the proof of this theorem in Section 3.1 (see also 
Remark 3.2), there are other choices for the coefficient c that yield critical vector¬ 


valued measures, but for which the associated Riemann surface is of genus 1. We 
do not consider this case here. 


Theorem 3.7 below (see also Remark 3.8) shows that the measures p a from 
Theorem 1 1.12 are in fact equilibrium measures on an appropriate family of contours 
in C. 


For the critical measure fi a given by Theorem 1.12 the algebraic function £ 


defined by ( 1.13| ) has two real branch points a± < b\ and two nonreal branch points 
& 2 , «2 = ^2 (Im 62 > 0). The support of the components /ij of the critical measure 
Pa = (/ii, /r 2 , p. 3 ) can be easily described; we show that it exhibits a phase transition 
for a certain value of a: 
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Figure 1. For T ss 0.126 < r c (left panel) and for r ss 0.227 > r c 
(right panel), numerical evaluation of the supports of the measures 
/ii (long dashed line), p .2 (continuous line) and [13 (short dashed 
line - only on the right panel). 


Theorem 1.14. There exists a critical value a c £ (0,1/2), determined by 


!m [ (£i(s) -€ 3 {s))ds = 0, 

Jbn 


(1.34) 


where ^,£3 are the functions in (1.12), such that 


• If a < a c then /Z 3 = 0, supppi = [ai, 6 i], and supp/i 2 is an analytic arc, 
disjoint from supp /ii and connecting the branch points 02 , b 2 . 

• If a c < a < 1/2 then none of the measures is zero, and there exists a value 
a» £ (ai,&i), determined by the condition 


Im / (£i(s) - £ 3 (s))ds = 0, 

Jb 2 


(1.35) 


for which supp /ii = [a*, 61 ], supp ^3 = [ai,a*] and supp /a* = 71 U 72 , where 
71 is an analytic arc on the upper half plane connecting 62 to a*, and 72 
its complex conjugate arc. 

We refer to Figure [T] for an actual numerical evaluation of the supports of the 
measures pi,ii 2 ,h 3 i n Theorem |1.14| 

Remark 1.15. The restriction a < 1/2 arises somewhat naturally in our analysis. 
For a = 1/2, the algebraic equation ( 1.13| ) with coefficients (1.321- (1.33) still defines 
a vector-valued measure /ti /2 in the sense of Theorem 1.3 but the support of its 
third component is unbounded; hence, strictly speaking the measure /ti /2 does not 
belong to M\/ 2 - 


In rough terms, Theorems |1.2| and |1.3| induce a map 

{critical vector-valued measures} 1 —> {algebraic equations} 

and Theorems mu and |1.14| characterize this map completely for a subclass of 
algebraic equations and the cubic potential (1.30). 

























16 


A. MARTINEZ-FINKELSHTEIN AND G. SILVA 


The existence of the critical measure p assured by Theorem 1.12 and the topo¬ 
logical description of its support given by Theorem |1.14 follow after a careful anal¬ 
ysis of the associated quadratic differential zu in (1.28). We describe the criti¬ 


cal graph of vo for a = 0 and indicate its subsets corresponding to the supports 
supp hi, supp p, 2 , supp/ 13 . In other words, we lift the supports to the associated 
Riemann surface, embedding them into the critical graph of vo. We then deform 
the parameter a in the interval ( 0 , 1 / 2 ), observing the dynamics of the critical graph 
and keeping track of the supports of the measures. 

During this deformation procedure, the critical graph of vo displays several tran¬ 
sitions, described in Sections [4.5.2 - 4.5.7 but only one of such transitions, governed 
by equation (1.34), has direct impact on the topology of supp/q, supp/q, supp/q. 


It should be stressed that there are not many tools for such an analysis (we have 
summarized the basic facts about quadratic differentials in the Appendix [B]), and 
that it can be extremely complicated, even on the Riemann sphere. The situation 
becomes even more involved when we consider trajectories on a more general com¬ 
pact Riemann surface. We believe that the methodology we have developed is of 
independent interest, and can be applied in other situations as well. 

1.4. Critical measures and max-min problems. The critical measure given by 
is crucial in the study of multiple orthogonal polynomial P n ,m(z) = 


1.12 


Theorem 

z N + ... defined through 

/ Z j P n 


l (z)e Nz3 dz = 0, j = 0,...,n-l, 


[ z J P n ^ m (z)e Nz3 dz = 0, j = 0, ...,m— 1, 

JY,n 


(1.36) 


where N = n + m and £i ,£2 are contours extending to oo with angles —^,0 
and — respectively. It turns out that the asymptotic behavior of P ram as 

m, n —> oo in such a way that n/N — )• a £ (0,1/2), is governed by this critical 
vector-valued measures p £ A4 Q : under this regime, the sequence of zero counting 
measures associated with (P ny m) converges to the sum of the first two components of 


p. The two extremal cases, a = 0 and a = 1/2, have been proved in 27 and 35 


respectively, and the general case a £ ( 0 , 1 / 2 ) will be addressed in a follow-up 


publication 59 


As we already mentioned, zeros of orthogonal polynomials on the real line are 
characterized as solutions to logarithmic equilibrium problems. The measure p 


given by Theorem 1.12 and hence the zeros of the multiple orthogonal polynomials 


(1.36), also admits a vector equilibrium problem characterization, as explained next. 


Let T be the class of contours £ extending to oo with angles — //, 2 ./ and inter¬ 
secting the real axis exactly once, say at a* = a*(£). Denote by AIo/E) the subset 
of measures v = (iq, zq, ^ 3 ) £ A4 Q that further satisfy 

suppzq c [a*, + 00 ), suppiq C £, suppzq C (— 00 ,a*]. 

The a-equilibrium measure v a = V a ,T. of £ £ T (for the energy functional E(-) 
given in ( 1 . 2 ) with ( 1.30|) (1.31)) is the unique measure v a £ A4 Q (£) for which 


E{V a ) = ^ inf E(y) =: E a {Y). 

OeM a ( s) 

Theorem 1.16. For a £ (0,1/2), there exists a contour T £ T for which its 
a-equilibrium measure is the a-critical measure p given by Theorem\l. 1 
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We conjecture that the contour T is characterized by the max-min property 

E a (T) = sup E a (E). 

£eT 

This characterization would be analogous to the max-min property for non- 
hermitian orthogonal polynomials [53p9f74|75] . Furthermore, an adaptation of the 
methods in 53 58 shows that the a-equilibrium measure of any contour solving 


this max-min problem is in fact an a-critical measure. 


1.5. Structure of the rest of the paper. Thcorems |1.2[|l .< 8 [| 1 ,9[ and a simplified 
version of Theorem 1.3 are proven in Section[2] (the complete proof of Theorem 1.3 
is rather technical and thus deferred to Appendix [A| . In Section 3.1 we derive the 


coefficients (1.32) (1.33) and prove the second part of Theorem 1.12 Assuming 


the existence of certain short trajectories (fact proved later, in Section |4j), we de¬ 
rive the first part of Theorem |1.12| and Theorem |3.7| in Section |3.2| which yields 
Theorem ra Sectio n [4] is devoted to the analysis of the global structure of the 
quadratic differen tial (|1.28[ ) for the cubic potentials (1.30), which also provides a 
proof of Theorem |1.14 (or equivalently, of Theorem 4.4 where the structure of the 
support of the critical vector-valued measure is rephrased in the terminology of 
quadratic differentials). Finally, two more appendices contain the minimal back¬ 
ground on quadratic differentials (AppendixjB]) and the description of the numerical 
experiments and procedures used in this paper (Appendix [C|). 

2. Critical vector-valued measures 
Critical measures were defined in the previous section in terms of the vanishing of 


the total energy (1.7). Here we use an alternative characterization, more convenient 


for calculations. When p = g is a scalar measure, the following Proposition 2.1 


coincides with 58 Lemma 3.1]. The proof extends trivially to vectors of measures 


p and polynomial external fields as considered here. We skip the details. 
Proposition 2.1. For any vector of measures p £ M a , 

E(p l ) = E{p) + t ReV h (p) + 0(t 2 ), 

where 

h(x) - h(y ) 


3 

T>h{p) = - a jt k 

j,k-l 


x-y 


-dnj(x)dn k (y)- 


3 

■E 

j=i ■ 


it is valid 
t -£ 0 , 

[ &Ax)h(x)dyLj(x). ( 2 . 1 ) 


Although the definition of critical measures only requires to evaluate T>h(p) for 
test functions h £ C 2 ( C), the quantity T>h(P) as defined above will be used for 
any function defined on supp pi U supp p 2 U supp /T 3 , as long as the integrals on the 


right-hand side of ( 2 . 1 ) exist. 

In particular, considering h and — ih , we get from Proposition ! 2.1 
Corollary 2.2. A measure p £ A4 a is critical if, and only if, 

V h (p) = 0 , 

for every function h £ C 2 (C). 


The following lemma is inspired by 58 , where a similar situation (but involving 


only one measure) was analyzed, and borrowed in the following form from 53 


Lemma 3.5], where a simplified proof is given. The proof from 53 is easily extended 
to our situation by simply choosing g = 1 therein and mimicking the arguments. 
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Lemma 2.3. Suppose jl is a critical vector-valued measure. If z € C satisfies 

/^<co, J = 1,2,3, 

J \x-z I 

t/ien 2?/^ (/2) = 0 for the meromorphic function 

1 


( 2 . 2 ) 


h z {x) = 


x G C. 


(2.3) 


The next result is straightforward, and we skip its proof: 


Lemma 2.4. Let f be liolomorphic in a punctured neighborhood of p £ C. If 
\f(z)\ r is mes 2 -locally integrable at p for some r > 0, then p is either a removable 
singularity or a pole of f of order < 2/r. 

For any finite Borel measure p, it follows from Tonelli’s theorem that the function 

/' df ^ 

J \x-z\ 


is locally mes 2 -integrable and in particular finite mes 2 -a.e.; that is, ( 2 . 2 ) is satisfied 
for mes 2 -a.e. 2 e C. Hence the Cauchy transform C M of the measure p is finite 
mes 2 -a.e. In particular, the functions £3 in (1.12) are finite mes 2 -a.e. and 

locally integrable with respect to mes 2 , so that 

£«?(*) (2-4) 

i=i 

is well-defined mes 2 -a.e. and is analytic on open subsets of C \ (supp p\ U supp p -2 U 
supp^ 3 ). 

Recall the quantities R, T>h{p) and h z defined in (1.14), and ( |2.3| ), respec¬ 
tively. The next lemma is crucial for what comes later: 


Lemma 2.5. If jl £ A4 a , then the functions (j defined in (1.12) satisfy the identity 


\ (£1 + £2 +£ 3 ) ( z ) = R i z ) +V hz (p) mes 2 -a.e. 


(2.5) 


Moreover, the statements 

(i) is a polynomial, 

(H) Ej=i^ 2 (-) = 2R ( Z ) mes 2 -a.e. 

(Hi) T>h z (fi) = 0 for every z satisfying (2.2), 
are equivalent. 

Proof. For h = h z , 

h z (x) - h z (y) 


x-y 


) ff 

-dpj(x)dp k {y) = - // 


(x - z)(y - z)“ n 
= -C«(z)C w («), 


dpj(x)dp k (y) 


and 


where 


J hzW&jitfdpjix) = &j(z)C^(z) - Qj{z ), 

QM-fW-W 

J X - z 


( 2 . 6 ) 


(2.7) 


-dpj(x). 


( 2 . 8 ) 
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Notice that the integrand in the right-hand side of ( |2.8| ) is a polynomial both in x 
and in z, and thus (since /i ? is finite and compactly supported) it is also dyj(x)~ 
integrable for all z. Hence, Qj is a bona fide polynomial. Additionally, (2.6) and 
(2.7) also show that Dh z (jl) is finite whenever (2.2) holds true. 

Using (2.6) and (2.7) in (2.1), we get that 

3 3 3 

ZU/Z) = E a iik C^{z)C^{z) + J2^) CH (z) - 

j,k =l i=i j=i 

or equivalently, 

3 

VhM = C(z) t AC(z) + ($'(z)) T C(z) - E Qj(z), (2-9) 

3 = 1 

with the previously used notation C = (C^ 1 , C M2 , C M3 ) T and i>' = ($' 1; <b' 2 , <b' 3 ) T . 

On the other hand, by Jl.lll) (and taking also into account both identities from 
(pOt) we get 


{& ++e 3 2 ) = l?? = +c^j a 

= l ($') t A$' + C T AC + % ($') T AC 

= l (d>') T A& + C T AC +\U') T (3/ 3 - bb T ) C 
y o \ / 

C t AC+ + *(4>') T A<1>' - i ^b T S'Y {b T <Y . 


Notice that the compatibility condition (1.3) means precisely that & T <N = 0, so 
that 

\ (£? + £ 2 + il) = l (&) t a$' +1 c T ac + (<£') T c|, 


and using ( |2.9| ) we can simplify the expression between brackets above to conclude 
that 

2 (£l + £2 + £3 ) = ®h z (fi) + g ( < 5 , ) T A& + E Qj = ’Dh z (m) + R, 

3— 1 

where R was defined in ( |1. 14 1. This identity is valid for values of z satisfying (2.2), 
which proves (2.5), as well as the equivalence of (ii) and (in). 

Clearly, (ii) implies (i). On the other hand, by (2.6), 

h z {x) - h z (y) 


lim 

Z—t OO 


x-y 


-d/jj(x)dfi k {y) = 0 , 


and since function h z (-) converges uniformly to 0 as z -A 00 , by the Dominated 
Convergence Theorem we conclude that 




lim / h z (x)^' i (x)dfxj(x) = 0. 

Z—t OO J J 

lim V hz ( jl ) = 0. 


Thus 


(2.10) 
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Measures yi, y -2 and y 3 are compactly supported, so that (2.5) holds for all z 
sufficiently large, and with ( | 2 . 10 [ ), 

lim (Ci+^ 2 +f 3 - 2 iZ)(«)= 0 . 


This establishes that (i) implies (ii). 


□ 


Proof of Theorem | ?.d[ A combination of Corollary |2.2| and Lemma |2.5| shows that 


if y is critical, then the polynomial R in (1.14) is alternatively expressed as 

Uzf). 


1 , 


R(z) = ^l(z ) 2 


■UzY 


From the expressions (1.12) it follows that 

£1 + £2 + £3 = 0 , 

and as a consequence 

R = 2 (£1 + £2 + £3) = £1 + £i£ 2 + £2 = £2 + £a£3 + £3 = £3 + £s£i + £?• 


( 2 . 11 ) 


( 2 . 12 ) 


Multiplying the identities in (2.12) by appropriate factors (£j — £*,), it follows 
that 


£i(A ) 3 - £1 (z)R(z) = £ 2 (A ) 3 - &{z)R(z) = £ 3 (z ) 3 - &(z)R{z). 


Define 


D (z ) = -^(z) 3 + ^(z)R(z). 


(2.13) 


(2.14) 


Due to (2.131, D does not depend on the choice of j £ {1, 2, 3}. When we choose 
j = 1, we see that D is analytic on C \ (suppp,i U supp/i 2 ). Similarly, when we 
choose j = 2, 3, we see that D should be analytic on C \ (supp y 1 U supp /i 3 ) and 
C \ (supp /r 2 U supp/i 3 ), respectively. Consequently, D can only have singularities 


at S a , defined in (1.5), so from ( 1 . 6 ) these singularities are all isolated and in a 


finite number. 

Moreover, 

i^)r /3 <(i£,-wi 3 +i£twn^)i) l/3 

<i£iWi+i£iwr /3 i^)i i/3 , 

where for the last inequality we used (x + y ) 1 ^ 3 < x 1 / 3 + y 1 / 3 , which is valid for 
any non negative numbers x, y. 

The function £j is locally mes 2 -integrable (see the comments after Lemma [2.4[ ). 
Hence the function D satisfies the conditions of Lemma 2.4 for r = 1/3 and any 
choice of p £ S a , so the points in S a are not essential singularities of D. Moreover, 
from the behavior of the £/s and R when z —» 00, it follows that D has polynomial 
growth at 00. We conclude that D is rational (although it could be identically 
zero). 

In summary, we have shown that under the assumptions of Theorem |1.2[ there 


exist a polynomial R, defined in (1.14), and a rational function D with possible 
poles in S a , such that the functions £ 1 , £ 2 , £3 satisfy 

£j(z) 3 — R{z)^j{z) + D(z) = 0, mes 2 — a.e. 

The fact that all the components /ij are supported on a finite union of analytic arcs, 
and that they are absolutely continuous with respect to the arclength measure of 


their supports is then a direct consequence of 19 Theorem 2]. 
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Let p be a pole of D. The functions £ 1 ,^ 2,£3 satisfy the cubic equation (1.13), 
so their local behavior near p is of the form 

ij{z) = {z - p)~ Vi {kj + °{ 1)), as z^p, j = 1,2,3, (2.15) 

where 1 ^ 1 , 1 ^ 2 , ^3 £ Q and k\,k, 2 ,K 3 are nonzero constants. Since p is a pole of D 


and R is a polynomial, from (2.14) we learn the u^-’s are all equal, say to v , and 


3u should be equal to the order of p as a pole of D. Moreover, we further get that 
none of the functions £ 1 , £2 and £3 is analytic near p, so p has to belong to the 
support of at least two of the measures /ii,/i 2 and /.13. On the other hand, since £ 7 
is a linear combination of polynomials and Cauchy transforms of finite measures, 
we must also have v < 1, so we conclude 


v £ < x 


1 2 
3’ 3 


The behavior (1.15) then follows from Plemelj’s formula and (2.15), concluding the 
proof. 

□ 


Remark 2.6. The fact that (scalar) finite measures, whose Cauchy transform sat¬ 
isfies a quadratic equation mes 2 -q.e., live on a finite set of analytic curves was 
without any additional constraint on the measure (the proof was 
This was extended to an arbitrary algebraic equation in 191, 


established in 
inspired by 


58 


14 


but with a crucial assumption on the measure that its support has a zero plane 
Lebesgue measure (this assumption is embedded in the definition of the class Ai a 
in Section 1.2). Nevertheless, we believe that the assertion about the support is 


still true even if we drop this a priori restriction on p. 


Our next goal is to prove Theorem |1.3[ In other words, for a vector of measures 
p £ M. a we assume that the components of p are supported on a finite union of 


analytic arcs, and the associated functions £ ? in (1.17) satisfy an algebraic equation 


of the form (1.13) for some polynomial R and some rational function D. As we 


already pointed out after the statement of Theorem 1.3 these assumptions allow 


us to write £ as in ( 1 . 12 ) for the vector of external fields = ($ 1 , $ 2 , $ 3 ) given by 


(1.18). In particular, it follows that £ 7 ’s satisfy (2.11) and consequently 


R = £16 + £16 + £26 = ^ (£1 + £2 + £3)- 


From Lemma 2.5 we get immediately the following result: 


Corollary 2.7. Under the assumptions of Theorem \l.S\ the vector of measures jl 
satisfies 

v h M = o 


for z £ C \ (supp pi U supp P 2 U supp ps), where h z is as in (2.31. 


Hence, th e pro of of Theorem [L3] is reduced basically to extending the conclusion 
of Corollary 2.7 to the whole class of functions h £ C 2 (C). Although this is rela¬ 


tively straightforward for h analytic in a neighborhood of the union of supp pj , it 
takes more effort beyond the analyticity due to the integrand ( h(x ) — h(y))/{x — y) 
in (2.1). In this case we need to appeal to a “refined” version of the celebrated 


Mergelyan’s theorem that would allow for an approximation of a given sufficiently 
smooth function by a sequence of analytic ones, and in such a way that the sequence 
of derivatives is uniformly bounded. The complete proof we were able to find is 
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rather technical, so for the sake of readability we postpone it to Appendix [A] and 
present here the arguments valid for the case when h vanishes in a neighborhood 
of the set of double poles of D (or when this set is empty). We formulate it as an 
independent proposition: 

Proposition 2.8. Let p £ A4 a be a vector of measures satisfying the conditions of 
Theorem 1.3 Suppose that h £ (7 2 (C) satisfies 

supp h D S a = 0, 

where S a is the set of double poles of D. Then Vh(P) = 0. 

Proof. Recall that the union of the support of the components of p is a finite union 
of analytic arcs, whose points of intersection are denoted by S a (see (1.5)), and the 
poles of the coefficient D in (1.13) belong to S a , as it follows from the arguments 
in (2.15) et seq. Furthermore, if p £ S a \ S a , then 

dp 


^ f(z) = 0(\z- P r ), as * 


■P, 3 = 1,2,3, 


(2.16) 


for some v < 1/3 (possibly v < 0). 
Given h as above, let us define 


g(z) = (gi(z),g 2 (z),g 3 (z)) T and H(z) = (Hi(z), H 2 (z), H 3 (z)) q 


with 


( ^ [ h(x) f h(x) - h(z) 

9j(z) = / —— - d Pj( x )y H j( z ) = / -—- dfj,j(x), 

J X Z J X z 


J = 1,2,3, 


where as usual we understand gj(z ) in terms of its principal value (as in (1.81), for 
z £ supp^j. With this convention and recalling (1.24|, on any open analytic arc 
£ from S Q , not containing any branch point of (1.13), by the Sokhotsky-Plemelj’s 
formula we have 


gj±(z) = ±nih(z)pj(z) + gj(z), zeE, 


(2.17) 


where pi = dpj / ds is the Radon-Nikodym derivative of pj with respect to the line 
element on £, so that if £ D supp Pj = 0, we have pi =0. 

Let £ be an open analytic arc from supp p^ U supp p 2 U supp p 3 , not containing 
any branch point of ( 1.13 ). Observe that the matrix B in ( 1.10 ) satisfies the identity 
— M 2 , where 


Mi = 



M 2 = 



Since H is continuous across £, by (1.10), 


z £ £. 


2AH(z) = (Mi - M 2 ) BH(z) = MiBH_(z ) - M 2 BH + (z), 

Furthermore, the identity 

H(z) = g(z ) - h(z)C(z) 

holds for all z £ C \ (supp pi U supp p 2 U supp p 3 ), and we obtain that 
2 AH(z) = (. MiBg_ - M 2 Bg+) (z) - h(z) ( MiBC _ - M 2 BC+ ) (z). 


(2.18) 


(2.19) 
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Recall that Equation (1.161 allows us to rewrite (1.17) as (1.12) with the external 


fields given in (1.18). From (1.11), 


BC = £--B<f>', 


so that 


MiBd-iz) - M 2 BC + (z ) = - ±B&) - M 2 (£ + - *£?<!') ) (z) 

= Mi|l( 2 ) - M 2 £ + (z) - 2 -A&{z). 


( 2 . 20 ) 


Because the functions £ 7 's satisfy a cubic equation, the boundary conditions 


(1.25) are valid and imply the equality 

*Mi£_(*) = xm 2 £ + (z), 


zeE, 


( 2 . 21 ) 


where 


X = X(z) = 


{z) 0 

0 X 2 (z) 

0 0 



and Xj is the characteristic function of supper,, j = 1, 2, 3. The matrix X is piece- 
wise constant; for ease of notation we suppress its ^-dependence in the following 
computations. 

Since 4>' is an eigenvector of A corresponding to the eigenvalue 3/2 (see Sec¬ 
tion 1.2), we learn from (2.20) and (2 .21| ) 

XM\BC-{z) - XM 2 BC+(z) = -X$'(z), zeE. (2.22) 

An immediate consequence of this identity, combined with Sokhotsky-Plemelj’s 
formula C + = +C-, is 

IMi BC+{z)-XM 2 BC-(z) = -X&{z)+2niX(M 1 +M 2 )Bfl'{z), zeE, (2.23) 

where as usual, /t = (g 1 , H 2 , H 3 ) T 1 etcetera. 

Multiplying (2.19) by X and using (2.22), we conclude that 

2 XAH(z) = X (MiBg^ - M 2 Bg+) ( 2 ) + h(z)X&(z), 2 G E. 


Since 

(M lB f = - M 2 B , X T = X, 
we can write it equivalently as 


(am) T x) (z) = (glM ± BX - g T _M 2 BX) (z) + h{z ) ($'(*)) T X, 2 


G E. 
(2.24) 

Recall that E C S Q was arbitrary, so we can integrate this formula along 
Slipp Hi U SUpp g, 2 U SUpp H3- 

By ( |2.17| ), for j e {1,2,3} we get 


9j±(y)dHk{y) = ±M / h(y)H'j(y)dHk(y) + / gj(y)dg k {y ), 


(2.25) 


where we recall that if gj is discontinuous across supp Hk, we understand the second 
integral in the right hand side above in terms of its principal value. By assumption, 
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supp h does not contain double poles of the coefficient D , and we use the condition 
(2.16) to get the behavior 

h(y)t i j(y) = °(\y-p\~ u )^ as y~*v ( 2 . 26 ) 

for some v = v(p) <1/3 and any endpoint p of the support of y k . In particular, 
this implies that /i/z' and gj are /r^-integrable, so the integrals in the right-hand side 
of (2.25) are convergent. Moreover, the behavior (2.26) allows us to interchange 
the order of integration for the second integral in the right-hand side of (2.25) 
(see 62 Equation (20), page 25]) in order to get 

J 9j{y)dpk{y) = - j h(y)C fJ,k (y)dfj,j(y). 

Again by Sokhotsky-Plemelj’s formula, 

f Ky)C llk {y)dn j (y) = zprri f h(y)y' k {y)dyj(x) + [ h(y)C£ k (y)dy j {y), 


and using it in (2.25), we finally conclude that 

J 9 j±(y)dy k (y) = ±2ni J n' j {y)lJ-' k (y)h(y)dy - J h(y)p, l j {y)C^ k (y)dy. 

A direct consequence of this formula is that for any 3x3 constant matrix Ad, 

J 9 ±Mdjl=± 2 TTi J h(y)(fl'(y)) T Mjl'(y)dy — J h(y)(jl'(y)) T MC±(y)dy. 

In particular, we have 

J ( g+MiB - g^M 2 B)dp, = 2t ri J h(y)(y'(y)) T (Mi + M 2 ) Bg\y)dy 

- J h{y){jl'{y)) T [M.BC+iy) - M 2 BC.{y )) dy 
= 2iri J h(y)(jl'{y)) T X(y) (Mi + M 2 ) Bjl'(y)dy 
- J h{y){it'{y)) T X{y) (M.BC+iy) - M 2 BG_(y )) dy 
= J Hy)(id'(y)) T ®'(y)dy, 

where we have used the trivial identity p! T X = fl ,T for the last two equalities, 
and also (2.23) for the last equality. Thus, integrating (2.24) with respect to dy, 
observing that Xdy = dy and applying this last identity, we arrive at 

J (ah) ( z )dy{z ) = J h{z) ($'( 2 )) dy{z), 

which is a compact form of writing the condition T>h(y) = 0. □ 


Proof of Theorem \ 1.8\ For a vector-valued function / we understand by 

J f(z)dz 

the term by term integration, and we denote by y’ = (y\, y’ 2 , y' :i ) the vector whose 
components are the densities of yi, y 2 , yz (with respect to the complex line element) 
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along the arcs of S a . Because the supports of the components of jl are made of ana¬ 
lytic arcs, and ft' has analytic components, we conclude that U = C7 M2 , U^ 3 ) T 
is continuous across the arcs of E Q . 

Integrating (1.11) and using (1.9) we get 

Re J £{y)dy = Re Q^B$(z) \ + BU (z ) - d, 

where d is a real constant vector which only depends on the connected component 
of C \ supp Hi U supp H 2 U supp H 3 - 

Let E be an open analytic arc from E a , not containing any branch point of 
(1.13). Reasoning as for (2.181, we have 

2AU(z) = ( M l - M 2 ) BU{z) = MiBU_(z) - M 2 BU+{z ), zeE. 

Thus, for z £ E and d± the vector of constants for the component on the ±-side of 
£, 

2 AU(z) = M 1 ^Re f £_( y)dy - Re + <L 


- M 2 ( Re / £+{y)dy - Re ( -B$(z) ) + d + 


= Re 


(Mi€-(y) - M 2 £ + (y)^ dy- - Re A$(z) + M x d- - M 2 d+. 


Recalling that A<^> , = (3/2)$' this last equation implies 

2AU(z) = Re J ~ M 2 £ + (y)j dy - Re$(z) + M x d_ - M 2 d + , z £ E. 

If E C supp Hi-ji then continuous across E, and (1.25) tells us that the (4— j)- 
th entry of Mi£_ — M 2 £ + vanishes, which in turn yields (1.26) (with j replaced by 

4 -j). _ 


Note that due to (1.26) the 5-property on supp Hi-j is equivalent to 

3 


d_ 

~dz 


Y a A- j , k U >lk 


\k =1 


^4-j 


d_ 

dz 


0‘A-i,kU llk + 


4>A-j 


\k =1 


(2.27) 


Equation (1.11) means that for z £ 


d ( 


1 


2- iAU(z) + -<f>(z) ) = AC(z) + -&(z) = -B T £(z), 


where we also used the identity involving A and B in (1.10) and A<&' = (3/2)<E>'. As 
before, the boundary conditions (1.25) imply that the entry (4 —j) of B T {£_ +£+) 
vanishes, and the equation above then implies (2.27). 

□ 


Proof of Theorem \1.9\ We denote the restriction of Q to the sheet IZj by Qj. In 
cyclic notation mod 3, Qj can be expressed as 

Qj(z) = (€j-i{z) - 0+i(^)) 2 ) z£U j , j = 1,2,3. 
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Clearly Q 2 is meromorphic in C \ (supp U supp p 2 U supp [ 13 ). We first prove 
that each Q 2 is meromorphic on the sheet IZj. If F C supp ^3, then due to (1.6) 
the function £1 is analytic across T, hence 

6±(s) = £3 =f00i s e r, 

implying that 

Q?+(«) = (6+00 - Ms)) 2 = (Ms) - Ms)) 2 = Ql-(s), s € r, 

so Q\ is analytic across F and, as a consequence, it follows that Q\ is meromorphic 
on the whole sheet 7Z\. Similarly we prove Q 2 is meromorphic on 7 Zj, j = 2,3. 

We now show that Q 2 is globally defined, that is, the function Q 2 +1 is the 
analytic continuation of Q 2 from 7 Zj to TZj+i, j = 1,2,3. 

Consider an arc T C S Q . For the construction of TZ we know that T connects 
exactly two sheets, which we assume to be 7Z\ and 72-2, the remaining cases are 
analogous. Then T C supp/iq \ (supp/1,2 U supp/x 3 ), so the function £3 is analytic 
across F, hence £i± = £2^ on F and 

Ql±(s) = (Ms) - M)) 2 = (Ms) - M)) 2 = QlM 

so Q\ is the analytic continuation of Q\ to 1Z 2 across T. 

Thus, Q 2 is meromorphic along any arc connecting the sheets. It is clear that 
the functions £i,£ 2,£3, and hence Q 2 , can only be unbounded at the points at 00 
and also at the poles of the coefficient D in (1.13). If p is such a pole, then the 
three functions £i ,£ 2 ,£3 are branched at p , and the local coordinate of 1Z at z = p 
is of the form z = p + u 3 , u £ C, and hence 

dz 2 = 9 u 4 du 2 . (2.28) 

Since p is at most a double pole of D, it follows from the identities 

D(z) = -Mf + R(z)U*)> J = i.2,3, 

that the functions £i,£ 2,£3 all blow up with the same order 1//3, v £ {1,2}, that is 

^( z ) = ^ 7 J^( l +°( 1 ))^ as z-+p, j = 1,2,3, 

for some nonzero constants k i, « 2, «3, so 
12 


Q Z ( Z ) = ^73 C 1 + °(!)) = + °( 1 )). 


as z 


for some constant c. Combining with (2.28), we get 

—Q(z) 2 dz 2 = cw 4 ~ 2l/ (l + o(l))du 2 , 

for some constant c. Since v < 2, we get that w has to remain bounded near p , 
and hence w can only have poles at the points at 00. 

Let E C 2 q FI supp pi be an open arc. Due to assumption (1.6), function £3 
is analytic across E. Now the definition of function £1 in (1.12), the Sokhotsky- 
Plemelj’s formula and (1.25) yield that for s £ E, 

dpi(s) = dpi(s) + dp 2 (s) = — (£i+(s) - £i-(s)) ds = — (£i+(s) - £2+00) ds. 

Analogously, if E C D supp ^2, then £2 is analytic across E, and for s £ E, 

dp 2 (s) = dpi(s) + dfj, 2 (s) = 


^ (6+00 - £i-(s)) ds = (6+00 - £ 3 +(s)) ds. 















CRITICAL MEASURES: STRUCTURE OF TRAJECTORIES 


27 


Finally, if S C S a fl supp/1,3, then 6 is analytic across E, and using again the 


definition of £3 in (1.12) and the Sokhotsky-Plemelj’s formula, we get that for 
s £ E, 

dii 3 (s) = dn 3 (s) - dm(s) = ( 6 + 6 ) - 6 - 0 )) ds = ( 6 + 0 ) - 6 + 0 )) ds. 

This shows that the measures fij must be supported on arcs of trajectories of w. □ 


3. The cubic case 

In the following two sections we deal with the cubic case, and describe a one- 


parameter family of critical vector-valued measures for the energy functional ( 1.11 
(1.2) and for the choice (1.30), so that the external fields are given by ( |1 .31 ). 
Although the cn-critical measures were defined for a £ [0,1], here we restrict our 
attention to a £ [ 0 , 1 / 2 ]. As it will follow from our analysis below, as a 6 1 1 / 2 , 
the support of the component p ,3 of the a-critical measure becomes unbounded. 
Furthermore, our original motivation was the asymptotic analysis of the multiple 
orthogonal polynomials (1.36), for which the case of a € [1/2,1] can be easily 


reduced to a £ [0,1/2] by an appropriate rotation of the plane. Thus, the selection 
of this interval for a is natural in the present situation, although our method carries 
over without any special difficulty to the whole range of a. 


3.1. The spectral curve. Recall that by Theorem |1.2| the shifted resolvents, de¬ 


fined in ( 1 . 12 ), of any critical vector-valued measure jl = (/zi, p. 2 , p. 3 ) satisfy the 


algebraic equation (“spectral curve”) (1.13). As a first step, we deduce the expres¬ 


sions (1.32)-(1.33) for its coefficients. 


For the potentials as in (1.31), the coefficient R in (1.14) reduces to 


R(z) = 3z 4 — 3z — c, (3.1) 

where at this moment the constant c is given in the form 

c=/ x(dfii(x) + dn 2 {x)). 

Since D(z) = —+ R(z )6 = 66(6 + 6 )> comparing the expansions of both 
expressions at 00 and using (|1.4|), we further get that 


c<“(*) + c»M = -j-2L 


1 — a(l — a) 
3 z* 




00 , 


and 


D(z) = — 2 z 6 + 32 s + cz 2 — 3r, r := a(l — a). 


(3.2) 


Observe that r = a(l — a) is an equivalent parametrization that gives a bijection 
between the interval a £ [0, |] and r G [0, \]. As it was mentioned in Section 

were studied in 


the extremal cases, r = 0 and r = 


27 


and 


in their connection to the multiple orthogonal polynomials (1.36iTfor the choices 


35 


1.4 


respectively, 


n = 0 and n = to, respectively, and the family of critical measures depending on a £ 
( 0 , 1 / 2 ) are also relevant to the asymptotic analysis of these polynomials for general 
n, m. From a different perspective, we are studying a continuous deformation of 
the critical measures, interpolating the extremal cases of r = 0 and r = /. 
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The discriminant of ( 1.13| ) with respect to the variable £ is 
Discr(z) = 4 R 3 (z) - 27 D 2 (z) 

= (81 - 324r)z 6 + 54 c 2 5 + 9c 2 z 4 + 54(9r - 2 )z 3 
+ 54c(3r - 2 )z 2 - 36c 2 z - 4c 3 - 243r 2 . 


Since Discr(z) is a polynomial in 2 of degree 6 , if the Riemann surface of (1.13) 
has genus 0 then, by the Riemann-Hurwitz theorem, Discr(x) must have a multiple 
root. In particular, the discriminant Discr! of Discr(z) must vanish. A cumbersome 
but straightforward calculation (that can be carried out with the aid of a symbolic 
algebra software such as Mathematica) shows that 

Discr i = arpi(c)p 2 (c ) 3 
for a nonzero real constant a and 

pi(c) = 64c 3 + 243(1 - 4r) 2 , 

p 2 (c) = c 6 - 486c 3 t(1 + r) + 2187r(3r - l) 3 , ’ ‘ ' 

so we expect c to be a root of either pi or p 2 . For the case r = 
the algebraic equation (1.13) reduces to 


j studied in 


35 


£ 3 - (3 z A - 3F)£ - 2 z 6 + 3 ~ 3 - - = 0 , 

showing that c = 0, which is a root of p\. By continuity, we expect c to be a root 
of p\ for every choice of r € ( 0 , j), concluding that 

3 243 2 

C =-64- (1 - 4T) ' 


or 


,' 243 ,i , \ 2 

: =-| 


(3.4) 


which is the same as (1.33). 


Obviously, c defined by (3.4) can take three possible values; for the rest of the 
paper we choose c to be real (and thus, negative), so that the algebraic equation 
(1.131 is also real. It should be pointed out that a different choice of c would lead to 
an algebraic equation corresponding to another triplet of critical measures jl, that 
can be obtained from the original one by rotation by ±27r/3. 

For the choice of c in (|3.4|), we can rewrite Discr(z) as 


943 

Discr(z) = — q 1 {z)q 2 {z) 2 , 


(3.5) 


where < 71 , g 2 are given by 


256 1/3 4 128 3 16 2/3 2 

- 4r ) 7 2 + -9-* + 3T73 (1 “ 4r ) 7 * 

- 3 1/3 32(1 - 4 t) 1/3 z + 16(1 - 8 r), 

q 2 {z) = 3 1 / 3 (l- 4 r) 1 / 3 2 -l. 


(3.6) 

(3.7) 


For r < v the discriminant of q\ i 


IS 


const xr(l — t) 2 (27 — 100r) 3 , 
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which never vanishes, and we conclude that q\ has always four distinct roots. For 
the choice r = |, q\ has two complex conjugate roots and two real roots, so this 
also holds for any r. 

The resultant of q\ and q 2 is 

const x (1 — 12t) 4 (1 — 4r) 2 ^ 3 , 

so that in the interval 0 < r < 1/4, the polynomials q\, q 2 share a root only when 
r = 1/12. It then follows that for r 7 ^ 1/12 the roots of <71 are branch points of 
multiplicity 2 of (1.131. 


If the double zero of Discr (that is, the zero of ( 72 ) is a branch point of (1.13), 
then its multiplicity (as a branch point) has to be three. This means that the three 
solutions to (1.131 should coincide for 2 = 6 *, so (1.13) has to share a root with 
its second ^-derivative, and hence fj(b*) = 0, j =, 1,2,3. Plugging this back into 
(1.13) we see that D(b*) = 0. But 

N (12t — l) 3 

^ ~36(1 — 4r ) 2 ’ 

so D(b*) = 0 only for r = 1/12. Hence for r 7 ^ 1/12 the point is a regular point 
of (3.9), that is, the algebraic equation (1.13) is not branched at z = 6 *. We already 
observed that the simple zeros of Discr are always branch points of multiplicity 2, 
so the Riemann-Hurwitz formula says that for r 7 ^ 1/12 the associated algebraic 
equation has genus 0. Continuity with respect to r assures that the genus cannot 
increase for r = 1 / 12 , and hence the genus is also zero for this value. 

The discussion above can be summarized in the following proposition: 

Proposition 3.1. The algebraic equation (1.13) with coefficients given by (|3.1|), 


(3.2) and (3.4), has four branchpoints (the zeros of the polynomial q\ in ( |3.6 | ) ) and 
a double point (the zero of q 2 in ( |3.7 | ) ). Two of the branch points are real and the 
other two form a complex conjugate pair. For r 1/12 all these points are distinct, 
while when r = 1/12, the double point and one of the real branch points of (1.13) 
coalesce. The associated Riemann surface has always genus 0. 

Although the choice of c in (3.4), as a root of p\ instead of p 2 , was mostly 
motivated by the construction of a continuous one-parameter family of critical 
measures /r interpolating the extremal cases a = 0 and a = 1/2 studied in 


27 35 


formulas (3.1), (3.2) and (3.4) can also be explained in terms of the genus 0 ansatz 
stated in Theorem 1.12| that we prove next. 

Proof of the second part of Theorem \ 1.1 S\ We already noted that if the Riemann 
surface has genus 0 , then c must be a root of at least one of the polynomials pi 
and P 2 in (3.3), and Proposition 3.1 assures that the roots of p\ give rise to genus 
0 Riemann surfaces. It thus suffices to show that zeros of P 2 give rise to a surface 
of genus 1 . 

Consider c to be a root of the polynomial p 2 - As before, c can be assumed to be 
real, the remaining nonreal choices of c as a root of p 2 can be reduced to the real 
ones with suitable change of variables. Hence, 

c = 3(9t+9t 2 + \/3(1+9t)t 1/2 ) 1/3 or c = 3(9t+9t 2 -V3(1+9t)t 1/2 ) 1/3 . (3.8) 
Making the change of variables r = u 6 / 3, this can be written as 

'3 \ 1/6 


c = c(u) = 3u(u + l)(ir — u + 1 ), 0 < |u| < 


V 
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where u > 0 corresponds to the first value of c in (|3.8|) and u < 0 corresponds to 


the second choice of c in (3.8). 


With this identification, the discriminant Discr(z) simplifies to 
Discr(z) = —27 [u + z) 2 q(z), 

where 

q(z) = z a (4u 6 - 3) - z 3 (8u 7 + 6m 4 ) + z 2 {9u 8 + 6m 5 ) 

- z(10u 9 + 12 m 6 - 4) + 5m 10 + 12m 7 + 12m 4 + 4m. 
The discriminant of q is given by 

— 6912(m 3 + l) 6 (10m 3 + 9) (3m 6 + l) 3 ^ 0 for |m| G (0, (3/4) 1/6 ). 
Moreover, 

q(—u) = 9m 4 (2m 3 + l) 2 , 

hence q{—u) 0, unless u = — 1 /2 1 / 3 , which corresponds to r = 1/12. But for this 


latter choice, the corresponding value of c in (3.8 1 coincides with the value of c in 


(3.4), hence the coefficients R and D are given by (1.32). 

Thus, let m ^ — 1/2 1 / 3 . Then Discr(z) has four simple roots, which have to be 
branch points of the equation, and one double root z = — u. At this double root, the 
algebraic equation ( |1.13| ) simply reduces to £ 3 = 0, so its three solutions coincide. 
This is only compatible with the fact that Discr(z) has a double root at 2 = — u if 
the three solutions are branched at this point. Hence we have four branch points 
of multiplicity 2 and one branch point of multiplicity 3, and the Riemann-Hurwitz 
formula gives us that the Riemann surface has genus 1. □ 


The first part of Theorem |1.12[ claiming existence of a-critical measures, will be 
given by Corollary 3.6 


Remark 3.2. Obviously the quadratic differential (1.28) still makes sense if c is a 


root of P 2 as in the proof above. Numerical experiments performed to compute its 
critical graph indicate that the corresponding algebraic equation should also give 
rise to a-critical measures as in Theorem [T3I 


3.2. Equilibrium problem from the spectral curve. In this section, starting 
from the algebraic equation 

e~R(z)£ + D(z) = 0, (3.9) 

where it is assumed that the coefficients R and D are given by ( |3.1[ )-( |T4| ), we find 
a vector of measures \x = H 3 ) € A4 a for which the respective ^-functions in 

(1.12) satisfy (3.9), and consequently we prove Theorems 1.14|and|1.16 A central 
object for this analysis is the associated quadratic differential (1.28) (see also (3.11) 
et seq. below) and its critical graph, whose description is postponed to Section |4j 
According to Proposition 3.1 the spectral curve (3.9) has two real branch points 


oq < bi, two non real branch points 62 = < 22 , with Re &2 > 0 , and a double point 

1 


6 * = 


(3(l-4r))!/3 


>0 for tG (0,1/4). 


gi(M = 


32(1 - 12t) 2 
9(1 - At) 


>0 for r G (0,1/4), 


Since 
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we easily conclude that 6 * > bi for r £ (0,1/4), with 6 * = bi if and only if r = 1/12. 

Moreover, equation (3.9) defines a three-sheeted Riemann surface 1Z of genus 
0 for every value of the parameter r £ (0,1/4). In contrast to (1.27), where the 
cuts for the Riemann surface are defined in terms of the supports of the critical 
measures, here the cuts that split TZ into the sheets TZi, IZ 2 , 72-3 can be chosen in a 
somewhat arbitrary way, as long as they connect the branch points. In our case, 
the sheet structure will be given by two cuts: the interval of the real line connecting 
the real branch points ai, b\ of 1Z and a Jordan arc (to be defined precisely below) 
that joins the complex-conjugate branch points 02 , 62 - This arc intersects R in a 
unique point a*, whose location depends on the value of the parameter r. Namely, 
there is a certain critical value r c ss 0.19, defined later, such that a* £ (ai,b\) if 
and only if r c < r < 1/4 (what we call the supercritical regime), see Figure [2] It is 
important to stress here that r c > 1 / 12 . 

Obviously, it ultimately follows from our arguments that this cut structure co¬ 


incides with (1.27) for the critical measure given by Theorem 1.14 

Using these cuts we define three oriented sets, Ai, A 2 and A 3 on C. Namely, 
A 2 is always the projection onto C of the cut joining the complex-conjugate branch 
points 02 , 62 and oriented from 02 to 62 - In the subcritical regime (0 < r < r c ) 
we denote by Ai = [ai,&i] C R oriented from ai to 61 , and set A 3 = 0. In the 
supercritical regime (r c < r < 1/4), Ai = [a*, b{\ C R and A 3 = [ai,a*] C R, 
both with the natural orientation. The orientation induces the left (denoted by the 
subscript “+”) and right (with the subscript “—”) boundary values of functions 
defined on C or K. 

We define 


7Zi — C \ (Ai U A 2 ), 72-2 — C \ (Ai U A 3 ), IZ 3 — C \ (A 2 U A 3 ), 


and build the Riemann surface TZ = TZ\ U 1 Z -2 U 7£ 3 , associated to the algebraic 


equation (3.9), connecting each pair of sheets 1Z 3 crosswise across the cuts as follows: 
Ai connects 1Z\ with TZ 2 , A 2 connects 7Z\ with TZ 3 , and A 3 connects TZ 2 with 7 Z 3 
(this last condition is clearly non-trivial only in the supercritical regime, when 
r c < r < 1/4), see again Figure [ 2 ] 

As we will see soon, the cuts Ai, A 2 , A 3 can be specified in such a way that we 
will be able to define positive measures living on Aj’s, and the construction of the 


sheet structure will coincide with (1.27). 


The three solutions £i,£ 2,£ 3 of (3.9) are enumerated according to their asymp¬ 
totic expansion at infinity: 


£ 1 ( 2 ) = 2 z 2 - 1 + 0 (z 2 ), 
z 

U*) = -** + -+0{z-% (3.10) 

Z 

&(z) = -2 2 + — +0(z- 2 ). 
z 


As before, the functions £ 1 , £ 2 , £3 are regarded as branches of the same meromor- 
phic function £ : TZ —> C, defined by (3.9). As it is rigorously given by Proposi- 
4.3 below, the function A is the restriction of £ to the sheet TZj of the Riemann 


tion 


surface 1Z. 
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Figure 2. Sheet structure in the precritical 0 < r < r c (left) and 
supercritical t c < r < 1/4 (right) regimes. 


Due to the explicit geometric description of the sets Aj’s and sheets IZj's, it is 
straightforward to check that with 

[ 6 (^)- 6 W> onUi, 

Q( z ) = l £i(z) - &(*)> on n 2 , (3.11) 

on n 3 , 

Q 2 extends as a single-valued meromorphic function on TZ, and 

w = —Q 2 (z)dz 2 (3-12) 


is the corresponding rational quadratic differential on TZ. The details are also carried 
out in the general setting of Theorem |1.9| in its proof. 

One of the main outcomes of the discussion in Section [d] (see Theorem 4.4 and 
Remark 4.11) is that we can choose the cut A 2 connecting a 2,&2 to coincide with 


the trajectory along which 

! J (£i+(s) - £, 3 +{s))ds = const, 


Re 


z £ A 2 


(3.13) 


Furthermore, A 2 can be extended to an analytic arc F from e 2 .”' 00 to e^’oo in 
such a way that 


and 

Im J (£i+(s) - £3 +(s))ds = const, z G T \ A 2 , 

(3.14) 

see Figure [3] 
Further, 

6 ( 2 ) -&(z) ± 0 , z £ T\ {a 2 , 6 2 }, 

(3.15) 

and 

Re J (£i+(s) - £ 2 +(s))ds = const, zeAi, 

(3.16) 


£i+00 - 6 +( 2 ) ¥= 0, zeA 1 \{a 1 ,() 1 }, 

(3.17) 
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Figure 3. Pictorial representation of the sets Ai, A 2 , A 3 and T 
in the supercritical regime. 


and in the supercritical regime t c < r < 1/4, also 


Re 


J ( 6 +(«) - ^ 3 +(s))ds = const, 2 G A 3 , 


(3.18) 


and 


(3.19) 

(3.20) 


6 + 00 -&+(s) ^ 0 , 2 6A 3 \{ai,o,}, 

for which we refer to Proposition |4.3| 

Finally, it also holds 

£1(2) - £2(2) > 0, 2 > 61, 

£ 2 ( 2 ) - £ 3 ( 2 ) > 0, 2 < min{ai,a»}, 

and in the precritical regime r < r c , 

£1(2) - £2(2) < 0 , 2 G(a*,ai), 

for which again we refer to Proposition |4.3| 

We define measures /ii, H 2 , M 3 on Ai, A 2 , A 3 , respectively, through the formulas 

d»i{s) = ^ (£i+00 “ £2 +(s))ds, s G Ai, 
dn 2 (s) = ^-(£i +00 “ £3 +(s))ds, s G A 2 , 

(3.21) 

if 0 < T < T c , 

. (£3 +00 - &+{s))ds, sGA 3 , if T c < T < 1/4, 

where ds denotes the complex line element on the respective arc. Due to their con¬ 
struction, the supports of the measures /u-i, yu- 2 , /J 3 satisfy the claims of Theorem |1.14| 
(see also Theorem 4.41. 
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Remark 3.3. Using the definition in (3.11) we can describe formulas (3.21) saying 
that we build pj from the values of Q on the sheet TZ^j, j = 1 , 2 , 3 . 


Proposition 3 . 4 . Expressions (3.21) define positive measures pi, p2, P3- 

Proof. We prove the statement for the supercritical regime r c < r < 1/4; the 
precritical regime, somewhat simpler, can be analyzed similarly. 

From (3.13) and (3.16) it follows that the measures pi, p 2 , 1*3 are real. Moreover, 
the densities of p±, p3 with respect to the complex line element ds are continuous 
and non vanishing in the interior of A l5 A 3 , so the respective measures do not 
change sign. As for /x 2 , its density with respect to ds is continuous when restricted 
to either the upper or the lower half plane, but not across R. 

We start now with p\. For x £ Ai, we can write 


Mi([zAD = 


(3.22) 


where 


^i(z) = [ (£1 (5) - £ 2 (s))ds, z eC\ ((-OO, bi] U A 2 ). 


By the asymptotic expansion (3.10), we have that 

ipi(z) = z 3 + O(z), z -A 00 , 


and since £1 — £ 2 does not change sign on [ 6 1; + 00 ), see (3.17), we conclude that 

ipi(x) > 0 , x > b\. 

The function ip\ vanishes at the point z = 61 with order 3/2. Since we already 
know that if 1 should map Ai to *R, from its order of vanishing we further get that 
the imaginary part of if i+(x) is negative when x £ Ai is sufficiently close to the 
endpoint b±. From (3.17), we conclude 

ipi+(x) £ *R_, x £ Ai. 


By (3.22 1 , this establishes the positivity of p\. 

Furthermore, the density dpL^/ds is the analytic continuation of dp.\/ds across 
A 2 , hence dp, 3 is also positive. 

Finally, in order to establish the sign of p, 2 , consider first x £ A 2 D <C_, C_ = 
{z £ C : Im(z) < 0}. Denote by A 2 [a 2 ,a;] the subarc of A 2 from a 2 to x. Then 

A* 2 (A 2 [a 2 ,x]) = (3.23) 


2m 


where 


i’2(z)=[ (£i+(s) - £ 3 +(s))ds, z £ C\ ((-oo,&i] U A 2 ). 

J ao 


From (3.14) and (3.15) we learn that Re( , i/ 2 ) is monotone on r(e 2 00 , a 2 ). The 
asymptotics 

ip2(z) = z 3 + O(z), 2 40 o, 

shows that 

Rei/> 2 (z) -£■ + 00 , z -£ 00 along T, 

hence Re(i/> 2 ) is strictly decreasing on r(“Tr 00 , a 2 ), = 0. Since ^2 vanishes 

with order 3/2 on a 2 , we conclude that 

ip2+(,x) £ zR + , x £ A 2 flC_. 













CRITICAL MEASURES: STRUCTURE OF TRAJECTORIES 


35 


Thus, by (3.23) measure p -2 is positive on A 2 DC_. 

The positivity of P 2 on A 2 D C + can be obtained by similar arguments or using 
the real symmetry of the density of /r 2 ; we leave the details to the interested reader. 
The proposition is proved. □ 

Our next goal is an expression for the Cauchy transforms of combinations of 
measures p 1 , /r 2 , /T 3 : 

Proposition 3.5. The Cauchy transform of the measures p±, /x 2 , p^ defined in 
(3.21) are related to the £- functions in (3.10) through 

C^{z)+C^{z) + 2z 2 =fi(z), zG C\(A 1 UA 2 ), (3.24) 

C^(z) + C^(z)+z 2 = -&(*), * € C \ (A x U A a ), (3.25) 

C^(z)-C^(z)+z 2 = -^(z), zG C\(A 2 UA 3 ). (3.26) 

In particular, the total masses of the measures p\, /i 2 , p 3 satisfy ( |1.4| ), namely 

M + M = 1, |Mi| + 1/^31 = a, |/r 2 | - |/x 3 | = 1 - a. 

Proof. This is a straightforward consequence of residues calculations. For instance, 

6+(s) - &-(«) 


c« W + c»w = T^ 


f 1+( .) + J_ f 

s - z 2m J A2 


-ds 


1 

2m 
= Res 


600 

s — z 

m. 

Vs — z 


ds 


Res 


m, 

\s — z 


= -2z 2 + £ 1 ( 2 ). 


where the closed integral is computed along a contour oriented clockwise which 
encircles Ai U A 2 and does not encircle z. 

Analogously, 


C ,tl (z) + C' ia (z) = - Res 


600 


s — z 


s = z — Res 


6 W 


, s = 00 


s — z 


= -Z 2 - 6 (A)> 

and taking the difference of these last two equations, 

c»’(z) - C>*{z) = -Z 2 + Cl (z) + Uz) = ~z 2 - Uz). 


Finally, (1.4) is a direct consequence of the asymptotic expansion (3.10). □ 


A combination of Proposition |3.5| with Theorem 1 1. 3| yields the first part of The¬ 
orem m namely 

Corollary 3.6. The vector of measures ft = (pi, p 2 , ps) G A4 a defined through 
(3.211 is a-critical for the potentials (1.31). 


Recalling (1.9), the potential U M of a compactly supported signed measure p for 
which C \ supp p is connected is the real part of a primitive of C^, that is 

U^(z) = ReJ C^(s)ds + c, zGC\swppp, (3.27) 
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where the constant is chosen so as to have 


lim ^Re J C ,1 (s)ds + = 0. 


Let us apply this to the measures Hi, /i 2 , ^ 3 given in (3.21). From (3.24), 

(\c^{s) + C^{s))ds+\z*-\b\= r^ds, z£ C\(A 1 UA 2 ), 

Jb ± 6 * Jb 1 

and from (3.25), 

+ C^(s))ds + ^z 3 -Ul = - [~ 6 (a)ds, z £ C \ (Ar U A 3 ). 

Jb 1 6 6 Jb 1 

Summarizing, 

[ (2C Ml (5) + C* 2 (s) + C M3 (s))ds + z*-b\ 

Jb i 

= f (£i(s) - 60))ds, 2 : G C \ (Ai U A 2 U A 3 ), 

Jb 1 

where we must use the same paths of integration in the left and in the right hand 
sides. By (3.27) we then conclude 

2U^{z) + U^{z) + U^(z) + ^{z)~h =Re [ ( 6 (a)- 6 (s))ds, z £ C\A, (3.28) 

Jb i 

for some constant Z 1; the external field <fi given in (1.31) and A = A!UA 2 UA 3 . 
Since the £j’s have purely imaginary periods, the right hand side above is well 
defined regardless the path chosen. 

In a completely analogous way, we get 

2U^ 2 (z) + U^ 1 (z) — U^ 3 (z) + <f>(z) — Z 2 = Re f ( 6 (a) - 6 (s))ds, z £ C \ A, 

J a 2 

2U^(z) + U^{ Z )-U»%z)-l 3 =Re f ( 6 (a) - 6 (a))da, z £ C \ A. 

J min{ai ,0* } 

_ _ (3.29) 

(3.30) 

(3.31) 


Combining (3.28) with (3.16), we conclude that 

2t/' Ltl (z) + /7 M2 (z) + IP 13 (z) + 0(z) — Zi = 0, z£ Ar. 


Moreover, from the first equation in (3.19) we also get 

2U fll (z) + U> 22 (z) + U fl3 (z) + (l){z)-l 1 >0, z £ (6i,+oo). 

Furthermore, in the precritical case r < r c and z £ [a*,ai), 

2U^(z) + U» 2 ( Z ) + U^( Z ) + ct>(z)-l 1 =Re [ "(6 + (s) - 6+(s))ds 

J bi 

+ Re f (6(s) -£ 2 (s))ds 

./ ai 

/>ai 

= / (6(«)-6(a))da>0, (3.32) 

where for the second equality we used (3.16) and for the final inequality we used 

(pOol). 
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Analogously, 

2U» 2 (z) + U^(z) - U^{z) + </>(z) -l 2 = 0, zgA 2 , 
2U^ 3 (z) + U^{z) - U» 2 (z) - fa = 0, zcA 3 , 

2U ll 3 (z) + U tll (z) — U fl 2 (z) — f 3 > 0, z < min{a»,ai}, 
2U^ 3 (a*) + C/ Ml (a*) - [A 12 (a*) = Z 3 . 


(3.33) 

(3.34) 

(3.35) 

(3.36) 


If t > t c , we evaluate (3.301,(3.33), (3.34) at the common point a* € AinA 2 nA 3 
and take differences in order to get 

fa — fa — fa = 0 . 


On the other hand, if r < r c , we use (3.32) to get 

2t/ Ml (a*) + C/ A ‘ 2 (a*) + IP 13 (a*) + 4>(a*) > fa. 


We now combine this inequality with (3.33) evaluated at z = a* £ A 2 and (3.36), 
and conclude 

fa > fa — fa- 

Hence, for both the precritical and supercritical cases we can define 

Z 3 := Zi — fa, 

and with this definition, it follows that 

2U^ 3 {z) + U^{z) — U' J ’ 2 (z) — fa > 0, z < min{a*,oi}, 

2 U» 3 (z) + U^(z) - U 112 {z) - fa = 0, zeA 3 . 

One more variational inequality is based on the fact that we deal with a critical 
trajectory of our quadratic differential. Recall that by (3.29), 


2 U» 2 (z) + U^(z) - U» 3 {z) + 4 ){z ) - fa = Reip 2 (z), z&C\ A, 


where 


V>2 (z)=[ (€i(s)-£ 3 (s))ds. 

J ao 


From (3.14), we know that i [> 2 is real-valued on T \ A 2 . Hence, 

2 U » 2 (z) + U ^ 1 (z) - U 113 (z) + </>(z) - fa =Mz), zgT\A 2 . 


Moreover, by (3.15), the real-valued function i/j 2 is monotone on each connected 
component of T \ A 2 . Analyzing at oo, we see that 

i/j 2 (z) = z 3 + 0(z), 

hence ip 2 tends to +oo along T. Since it is zero at the endpoints a 2 , b - 2 of the 
connected components of V \ A 2 , we conclude that 4 > 2 is always positive on F \ A 2 , 
and hence 


2U> l 2 {z) + U^{z) - U^ 3 (z) + </>(z) -fa>0, z£T\A 2 . 


We summarize our findings in the following theorem: 


Theorem 3.7. Let measures fi\, /i 2 , /z 3 be defined in (3.21). 
constants fa, fa and 


Then there exist real 


fa fa — fai 
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h = 0, 

zeAi, 

V 

0 

2 e Ti \ A 3 

h — 0, 

z e a 2 , 

h > 0, 

z g r 2 \ a 2 


(3.37) 


such that the following variational conditions are satisfied: 

2U 111 (z) + U^ (z) + U* 3 (z) + </>(z) - h = 0, 

21/^ 1 (z) + U^ (z) + U* 3 (z) + </>(z) - 
2U» 2 (z) + U^ {z) - U» 3 {z) + <f>(z) - 
2t/ M2 (z) + U(z) - U» 3 (z) + </>(z) - 
2U fl3 (z) + (z) - U M2 (z) -l 3 = 0, zeA 3) 

2W* (z) + U^ (.2 ) - U 112 {z) -l 3 > 0, 2 g r 3 \ A 3 . 

where r 3 = [a*, +oo), T 2 = T, r 3 = (—oo,a*]. 

Remark 3.8. In the pre-critical case (r < r c ), A 3 = 0, thus the equality on A 3 in 
(3.371 is void. 

The equalities in ( 3.37| ) are the same as those in (1.26 1 . The extra information in 
Theorem |3.7| is coming from the remaining equations, which assure that the triplet 


(pi, p. 2 , M 3 ) is th e (unique) minimizer of the energy functional (1.2) over measures 
(^ 1 ,^ 2 , v 3 ) satisfying (1.4) (with mj replaced by zy,-) and supp/iy C Tj, j = 1,2,3, 
see [13] Theorem 1.8]. This is equivalent to Theorem 1.16 with T = T 2 . 

4. Global structure of the trajectories in the cubic case 

4.1. Dynamics of the singularities. A natural first step in the study of the 
structure of the trajectories of a quadratic differential is to clarify the position and 
the character of its singular points. In the case of the quadratic differential (3.12) 
we have to analyze the location and the dynamics of the branch points and the 
double point of the Riemann surface 1Z (corresponding to the spectral curve (3.9)) 
as functions of the parameter r. 

Recall that in Section 3.2 we have denoted the branch points of (3.9) by aq, a 2 , b\ 
and 62 , with the conventions 

ai, 6 iGK, ai < 61 , 02 = 62 , Ima 2 < 0 , 


and the double point of (3.9) by 6 *. 

Proposition 4.1. The main parameters of the Riemann surface TZ corresponding 
to the algebraic equation (3.9) exhibit the following behavior: as r grows from 0 to 

1/4, 

(i) the coefficient c in (3.4) increases monotonically from — 3 5 / 3 /4 to 0; 

(ii) • the branch point a\ decreases monotonically from 3 2 ^ 3 /4 to — 00 ; 

• the branch point b\ increases monotonically from 3 2 / 3 /4 to 3 2 / 3 /2; 

• the double point increases monotonically from 3 ~ 1//3 to + 00 . 

(iii) Always 

b\ <6*, 

and the equality is attained only for t = 1/12. 


Proof. Recall that we chose c in (3.4) to be real, thus (i) follows directly from the 
definition of c. 

Any branch point z = ( of 7Z is a zero of the polynomial <71 of degree 4 and real 

d T qi(0 


coefficients, defined in (3.6), so that 

drC = - 


d z qi(0' 


(4.1) 
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It is easy to check that the resultant of q\ and d T qi (with respect to the variable 
z) is 

const x (1 — 4t)“ 3(27 — 100r ) 3 ^ 0 for r G (0,1/4), 

which implies that for any root f of qi, <? t 9i(C) 7^ 0, and thus preserves its sign in 
the whole range of values of r £ (0,1/4). 

On the other hand, since ai (resp., b\ ) is the smallest (resp., largest) real root 
of q\, and the leading coefficient of q\ is positive, we know that 

d z qi(a i) < 0, d z q 1 (b 1 )>0. 


Straightforward calculations show that for r = 1/8, a± = 0 < bi and 

/ o \ 2 / 3 

d T q\{z) = —256 u 4 — 64 ((u — l ) 2 + l) , u = ( - j z. 

In particular, d T qi{z) < —64, and we conclude that in this case, 
d T qi(a i) = -128, d T qi(bi) < 0 . 


Hence, d T qi(a\), d T qi(bi) < 0 for all r G (0,1/4), and by (4.1l, a\ is a decreasing 
and b\ is an increasing function of r. 


From the expression (3.71 it follows that 

5 * 


! 


3 1 / 3 (1 — 4t)!/3 ' 


Replacing it in (3.6) we get 


qi{b *) = > 0 for t G (0,1/4), t ^ 1/12. 


9(1 - At) 


Since for r = 0, 


3 2 / 3 


ai 


= hi = — < 6 * = 3- 1/3 , 


(4.2) 


this concludes the proof of (iii). Finally, it remains to observe that for r = 1/4, 

, , 128 / 3 9\ 

which has one real positive root (3 2 // 3 /2) and two complex conjugate ones. This 
shows that 

3 2 /3 

lim ai = — oo, lim b\ = ——. 
r—>-1/4— r—fl/4— 2 

□ 


We finish this section with a technical lemma that will be used later. 


Proposition 4.2. For 0 < r < 1/12, the polynomial D in (3.2) does not have 
zeros on (ai,f>i), whereas for 1/12 < r < 1/4, D has exactly one zero on (ai,&i). 
Moreover, for 0 < r < 1/12, D(b\) = 0 only for t = 1/12, and D(ai) is never zero. 


Proof. On one hand, the discriminant of the polynomial D with respect to z is 

/(t) = const t( 1 — 12r ) 2 (l28r 2 — 32r — l) 

which shows that for r G (0,1/4), polynomial D has no multiple roots as long as 

1 / 12 . 










40 


A. MARTINEZ-FINKELSHTEIN AND G. SILVA 


On the other hand, we have seen in Section [3d] that the branch points of 1Z are 
the zeros of the polynomial q\ defined in (3.6 1 . The resultant (also w.r.t. z) of D 
and q\ is 

g{r) = const (1 — 12 t) 3 (27 — 100r) 3 , 


so again, for r £ (0,1/4), polynomials D and q\ have no common roots as long as 
r ^ 1/12, in particular implying that D does not vanish at oi, b± for r ^ 1/12. For 
r = 1 / 12 , we compute b\ = 6 * = 2 -1 / 3 and factor 


D(z) = -(z - b*) 2 (2z 4 + 2 5/3 z 3 + 3 2 1/3 z 2 + 2 + 2- 4/3 ), 


so for r = 1/12 we have D{b\) = 0 and D(a\) < 0. 

It is worth pointing out that both / and g can be easily found by means of any 
computer algebra software. 

Having in mind that D and q\ do not share roots for r ^ 1/12, it is sufficient to 
establish the assertion concerning the zeros of D for a single value of r in (0,1/12), 
and for a single value of r in (1/12,1/4). 

For r = 1/4, D is a quadratic polynomial in z 3 so its roots can be explicitly 
computed; we get that its only real roots are 


Zl 


M 


1/3 


Z2 


M 1/3 


Since 




from Proposition 4.1 (ii) we see that for r 

e > 0 , 


1/4 — e, for a certain small value of 


a\ < zi < bi < Z2- 


On the other hand, it is easy to check that for r = 0, 
D{ ai ) = D(b i) = D{ 3 2 / 3 /4) < 0, 


so that for r = e, for a certain small value of e > 0, D does not vanish on (ai, 6 i). 

□ 


4.2. The Riemann surface associated to the algebraic equation. In Sec¬ 
tion 13.21 we described the construction of the three-sheeted Riemann surface 7 Z = 
TZi U 7^2 U 7^-3 and of the branch cuts in such a way that the three solutions of 
(3.9), specified by the asymptotic conditions (3.10), become meromorphic on the 
respective sheet 7 Zj, with poles only at z = oo. They are also pairwise distinct as 
long as Discr(z), defined in (3.51, does not vanish, i.e. for 2 </ { 01 , 61 , 02 , 62 , 5 *}. 

The arc A 2 intersects K in a unique point a*; there is a critical value r c £ 
(1/12,1/4), to be specified later, such that a* < ai for 0 < r < r c (what we called 
the precritical regime), and a* £ (ai, 61 ) for r c < r < 1/4 (the supercritical regime). 

As a first result we establish some relations between the solutions £j(z). 


Proposition 4.3. Let t £ (0,1/4), r 7 ^ r c . Then 
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(i) for x £ R. \ (Ai UA 2 U A 3 ), 

66 ) < 66 ) < 66 ), x < min(a*,ai), (4.3) 

66 ) < 66 ) < 66 ), £>&*, (4.4) 

66) < 66) < 66), a* < x < ai, if t < t c , (4.5) 

66 ) < 66 ) < 66 ), bi<x<K, if 0 < r < 1/12, (4.6) 

£ 3 ( 2 ) < 66 ) < 66 ), bi < x < b*, if 1/12 < t < 1/4. (4.7) 


Additionally, 

6(6*) < 6(6) = 6(6), /° r 0 < r < 1/12, 

6 ( 6 *) = 6 ( 6 *) < 6 ( 6 *), /or 1/12 < r < 1/4. 
(ii) on Ai U A 2 U A 3 , 

o 

• for x £ Ai := A 1 \ {max(ai, a*), 61 }, 


660 = 606 e c\m, 66) sK, 

(4.9) 

and 


6± 6) = 6+ 6), 6+ 6) = 6- 6) ■ 

(4.10) 

0 

for z £ A 2 := A 2 \ {a 2 , 6 , max(ai, a*)}, 


6±6) = 6+6), 6+6) = 6-6)- 

(4.11) 

0 

for x £ A 3 := A 3 \ {a±, a*} (when r > t c ), 


66) =66) £ c\m, 66) e®, 

(4.12) 

and 


6±6) = 6+6), 6+6) = 6-6)- 

(4.13) 

Moreover, 


661) =661), if x < T c , 

(4.14) 

661) =661), if x > T c , 

(4.15) 

6(61) = 6(61), 

(4.16) 

662) =662) and 6(62) =6(62)- 

(4.17) 


Proof. Recall that we deal with the case a £ (0,1/2), so that 0<a<l — a<l. 
The behavior at infinity i n (|3. 10 > gives (4.3) and also (4.4). Furthermore, we have 
established in Proposition |4.1| that if r ^ 1/12 then &* ^ {ai, 61 , <3,2, 62 }, and in this 
case 6 * is a double zero of Discr(z). This means that when r yf 1/12, only two of 
the three values 6 ( 6 *), 6 ( 6 *), 6 ( 6 *) coincide, and since 

6 ( 6 *) + 6 ( 6 *) + 6 ( 6 *) = 0 , 


the coincident two differ in sign from the third one. Furthermore, since also 


( 666 ) ( 6 *) = -D( 6 *), 

we see that the sign of this third one is opposite to the sign of D(b *). But direct 
calculations show that 


D{ 6 *) 


(12r - l) 3 



for 0 < t < 1/12, 
for 1/12 < r < 1/4, 


36(1 - 4t ) 2 
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which means that for 0 < r < 1/12 the two coincident values of £j( 6 *) are the 
largest two, and the other way around if 1/12 < r < 1/4. It remains to use (4.4) 
in order to establish (|4.8[ ). 

Inequalities (4.6 H|4~7| follow by noticing that 6 * is a double zero of Discr(z 


so 


it is a simple zero of £i — £3 (for 0 < r < 1/12) and of £2 — £3 (for 1/12 < r < 1/4). 

From Proposition 4.2 we know that D(ai) ^ 0, D(b\) 7 ^ 0 for r 7 ^ 1/12. In 
particular, since for r = 0 (see (4.2)), a\ = b\ = 3 2 / 3 /4 and D(a\) = D(b\) = 


—81/2048 < 0, we conclude that D(ai) < 0 and D(b±) < 0 for 0 < r < 1/12. In 
other words, the smallest two of the three solutions £j come together at ai and at 
&!■ Now (4.6) implies (4.16) for 0 < r < 1/12. 

D(a 1 ) 


—00 as r 


1/4, 


On the other hand, recall that by Proposition |4.1 
and that for r = 1/4, 

D(z) = —2z 6 + 3z 3 — bi = 3 2 / 3 /2, 

so that D(bi) = 3/32 > 0. Using the same arguments we conclude that for 1/12 < 
r < 1/4, the largest two of the three solutions £j come together at b\ (and (4.7) 
yields (4.16)), and the smallest two become confluent at ai (and (4.151 follows from 

©). 

The discriminant Discr(:r) is negative for ai < x < b±, so just one of the solutions 
£j is real on aj < x < b 1 , and the other two are complex-conjugates of each other. 


Since we have ruled out the coincidence of the three branches, equality (4.16 1 implies 


that £ 1,^2 are non-real on Ai, which yields (4.9)—(4.10). By the same argument we 


also have (4.14), as well as (now using (4.3)) the identities (4.12) (4.13). 

Recall that we already established that the smallest two of the three solutions 
£j come together at ai, and (4.14) shows that for r < r c , 

£i(z) < £30*0. £ 2 ( 2 ) < £30*0, for a. <*< 01 . (4.18) 

For t < r c , when crossing A 2 two of the three solutions £ 7 are swapped, and the 


third one remains invariant. But the only option compatible both with (4.3) and 


(4.18) is (4.5). 


Finally, for r < r c , the inequalities (4.3) and (4.5) show that £2 is continuous 
across A 2 , and as a consequence (4.11) has to hold, obviously for the full range of 
t € (0,1/4). This, in turn, implies the last identity (4.17). 


□ 


Proposition 4.1 gives the formal proof of the fact that the Riemann surface 1Z, 


described in Section 3.2 is actually the Riemann surface of the cubic equation (3.9), 


and with this construction the function £j is meromorphic on 7 Zj and satisfies the 


asymptotic expansion (3.10). Moreover, £ : TZ —» C, given by £ = £y on 7 Zj, 


is meromorphic on 7 Z, giving the global solution to the algebraic equation (3.9). 


From Theorem 1.12 or alternatively Proposition 3.1 the Riemann surface 1Z has 
genus 0 . 

Let 7 r : 1Z —^ U be the canonical projection on the Riemann surface 7 Z. For a 
point p e C \ (Ai UA 2 U A 3 ), we denote by p^\ j = 1,2, 3, its preimage by n on 


IZj, that is, 


{p(U} = 7T 1 ({ p })nTZ j , j = 1,2,3. 

This notation is trivially extended to p £ Ai U A 2 U A 3 by taking boundary values. 
Notice that for the branch points it is valid 


b (1) ~ b (2) 

°i ~ u i 1 


b^=b^ 
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and 


a‘ 1 >=a‘ 2 > 


T> - -<» (r > r«), 


‘1 — “1 ( T A 7c), 

whereas if p belongs to a cut connecting exactly two sheets, say TZj and 7 Zk, then 

(j) ( k) 

p± =p t ■ 

We insist that the construction of 7Z is independent of the concrete choice of 
the cut A 2 ; this freedom will be used latter to specify an appropriate A 2 . Namely, 
in the next section we will show that A 2 can be made coincident with a critical 


trajectory of the quadratic differential w, defined in (3.12 1 ,which connects 02 and 
& 2 , see Definition 14.51 below. 


4.3. Computation of width parameters. Certain integrals of the function Q 


defined in (3.11) will play a crucial role in the upcoming analysis of the dynamics 
of the trajectories of the quadratic differential zu = — Q 2 (z) dz 2 on 7Z. They can be 
also formulated in terms of certain Abelian integrals on 1Z. 

Namely, we are interested in 


oj\ = u>i(t) = Re 


w 2 = w 2 (r) = Re 




U) 

,( 2 ) 


Q(s)ds = Re / (£ 2 (s) - &(s))ds, 


I b 2 


h( 2 ) 


Q(s)ds = Re / (£i(s) - £ 3 (s))ds, 


lb 2 


.(3) 


(4.19) 


w 3 = w 3 (r) = Re / Q(s)ds = Re / (£i(s) - &(s))ds, 


U3) 




0 J 4 = W 4 (t) = Re 


r-b* 




Q(s)ds = Re / (^ 2 (s) - £ 3 (s))ds, 


'b 2 


with w 4 (t) defined for r > 1/12. In this definition we understand that we integrate 
between two points on a sheet 7 Zj along a path that stays entirely in IZj. 

The values ojj are correctly defined regardless of the precise choice of the inte¬ 
gration paths. Indeed, the residues of the functions £i,£ 2 i £ 3 a t 00 are rea l (and 


independent of the value of r), see (3.10), so the integral of Q along a big loop on 
either IZj encircling 00 ^ is purely imaginary. It remains to notice that the genus of 
7 Z is zero, so that any closed contour around either branch cut A j can be deformed 
to such a big loop. 

Analytic computation of Wj’s is a formidable task. Instead, we have computed 
them numerically, see Figure [4] for the result. Since the integrands in ( |4.19| ) are 
multivalued functions, the numerical integration requires to implement an analytic 
continuation of each branch of £. We give further details in the Appendix [C| 

The functions wi, ui 2 , have two and one zeros on (0,1/4), respectively, whereas 
w 3 and u >4 do not vanish on the intervals (0,1/4) and (1/12,1/4), respectively. 

We denote the zeros of wi by n, t 2 and of 0 J 2 by r c . We have 


ri « 0.12487351, 


0.1913565, r 2 sa 0.2289555, 


so that they satisfy 

0 < ^ < n < r c < r 2 < 1/4. (4.20) 

We point out that the value t c is the one used in Section [3~2] in the construction of 


the Riemann surface 7Z , and it is equivalently determined by (1.34). 
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Figure 4. The graphics of the functions wi (upper dashed line), 
u >2 (short dashed line), W3 (dashed line with dots) and 0 J 4 (contin¬ 
uous line). 


It is also clear from Figure [4] that uji 7^ 0 J 4 for r > 1/12. 


4.4. Critical points of the quadratic differential. We follow the construction 
) and define the meromorphic quadratic differential w = —Q 2 (z)dz 2 
on the Riemann surface TZ. 

We start by analyzing the character of the critical points of w. For instance, the 
local parameter at z = is 


(1.28), (3.12 


z = a 


(i) 


+ « 


so that dz 2 = Au 2 du 2 . In consequence, {S,j—£,k) 2 {z)dz 2 has a double zero at z = 
if and only if £j( a i 1 ' ) ) 7^ ^(a^). A similar analysis at the rest of the points of 1Z 
yields the following classification of the critical points of vo on the Riemann surface 
(see Figure [5]): 

(a) For 0 < r < t c : 

(i) Simple zeros at a^\ b^\ a^\ ; 

(ii) Double zeros at b^ (only for r 7^ 1/12), a^\ b^; 

(iii) Double zero at b { 2 ' 1 if r < 1/12 and a double zero at b^ if r > 1/12; 

(iv) Zero of order 4 at = bi 1 ' 1 if r = 1/12; 

(v) Double pole at oc/ 1 ) with a real residue; 

(vi) Poles of order 8 at oc/ 2 ), 00 

(b) For t c < t < 1/4: 

(i) Simple zeros at a!{ \ b^\ ab^; 

(ii) Double zeros at a[ \ b^\ a£\ 6^; 

(iii) Double zero at bi 1 ' 1 ; 

(iv) Double pole at oc/ 1 - 1 with a real residue; 

(v) Poles of order 8 at 00 ( 2 \ 00 
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Figure 5. The finite critical points of vo in each case. Simple and 
double zeros are represented, respectively, with squares and dots. 


It is instructive to think of the critical points as evolving dynamically with r. 


Under this perspective, Propositions 4.1 and 4.3 show that for r small the double 


point (node) 6* corresponds to the double zero b^ on the second sheet. When 
t = 1/12, the points &*, b\ coalesce, giving rise to a higher order zero at the branch 
point = b^. For larger values of r, the double point emerges on the first sheet: 
now b l 2 '* is a regular point whereas &1 1 '* is a double zero. 


In the same spirit, the simple zero aq of the discriminant of (3.9) carries two 


(3) 

critical points of vo. For r < r c these are the simple zero a) ' and the double zero 
a'l' 1 = a{ \ For values of r larger than r c , these points interchange their roles: 
is a simple zero and = a\ 2> is a double zero. 

4.5. Analyzing the global structure of trajectories. 

4.5.1. General principles. The rest of this section is devoted to the description of 


the critical graph of the quadratic differential (3.12) for the whole range 0 < r < 


1/4. One of the outcomes of our analysis is the following theorem: 
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Theorem 4.4. For the quadratic differential zo = — Q 2 {z) dz 2 and for all values of 

( 2 ) ( 2 ) 

the parameter 0 < r < r c there exists a critical trajectory of vo joining a 2 and b 2 
on IZ 2 , whose projection by n on C is a real-symmetric analytic arc A 2 joining 02 


and 62 . 


( 2 ^ 

For t c < t < 1/4, there is an arc of critical trajectory of vo joining a 2 ' with a 
point a\ ; on the interval [a) ,b\ ; ] which is determined by (1.35), and the conjugate 
symmetric arc of trajectory joining b^ with the same point. The projection by tt 
on C of the union of these two arcs of trajectories is also denoted by A 2 . 


Recall that r c was formally introduced in Section [4.3| In virtue of the results in 
Section [3~2j in particular Corollary |3.6[ Theorem |4.4| implies Theorem |1.14| 

We remind the reader that up to now the branch cut, separating the sheets TZi 
and 7 ^. 3 , was free (see Figure [ 2 ]). In what follows we agree in the following: 


Definition 4.5. The curve, connecting 02 and 62 as part of the branch cut sepa¬ 
rating the sheets IZi and 72-3, is always given by the lift of A 2 (from Theorem 4.4) 
to the sheets TZi and 1Z 3 . 


In the next sections we will show that this definition is consistent with our 
construction of the Riemann surface TZ. 

One important fact is that the residues of zu at the poles at infinity (and the 
local behavior of the trajectories there) are independent of r: at oo^ they are 
closed analytic curves (so that 00 ^ is the center of a circle domain, see Appendix 
|b| , while oo( 2 ),oo( 3 ) attract trajectories in 6 asymptotic directions, given by the 
angles 

e loo ) = 2l-l^ j = 1 , e. (4.21) 


Critical values (4.20) split the interval (0,1/4) into the subintervals (0,1/12), 
(1/12,n), (ti,t c ), (t c , T2) and (t 2 ,1/4). We will show that the topology of the 
critical graph remains invariant in each of these intervals. 

The methodology we use can be summarized as follows: 


(i) Compute the critical graph for r equal to one of the critical values (4.20). 

(ii) Analyze the possible deformation of the trajectories for the values r + e, with 
e > 0 small, identifying the trajectories that display a phase transition. 

(iii) Prove that the topology of the critical graph is invariant inside the subinterval 
of interest, by analyzing the behavior of the widths wf s and showing that the 
corresponding strip and ring domains can not disappear. 

Along the way, will use some general principles that we enumerate here: 


P.l Quadratic differential vo has no recurrent trajectory for any value of r, see 
Jenkin’s Three Poles Theorem in Section fB~3l 
P.2 If 7 is an arc of trajectory of zu, then 7 , corresponding to the lift of the 
complex conjugate of 71(7) to the same sheet, is also an arc of trajectory. 
P.3 The complement of the critical graph of vj on TZ cannot have a simply 
connected component without poles on its boundary: that would contra¬ 
dict Corollary |B.3| or the maximum principle for harmonic functions on a 
compact Riemann surface. 









CRITICAL MEASURES: STRUCTURE OF TRAJECTORIES 


47 


P.4 The meromorphic function Q 2 * depends analytically on the parameter r. 
Hence the critical graph of w 1 and in particular all its critical trajecto¬ 
ries, depend continuously (in any reasonable topology, for instance, in the 
Hausdorff distance) from r. 

P.5 If for a certain value r = A, the point p belongs to the half plane domain 
for determined by the angles then the same holds true for 

a small neighborhood of values r £ (A — e, A + e), e > 0. The point p is 
also allowed to depend continuously on r. 

P.6 If for a certain value t = A, an arc of trajectory emanating from a given 
point p intersects the real line at a regular point, then the same holds true 
for r £ {A — e, A + e), e > 0. As before, the point p is allowed to depend 
continuously on r. 

There will be one more useful tool that we will employ several times in our 
analysis, formulated as Proposition |4.7| below. 

When describing the structure and the evolution of the trajectories of the qua¬ 
dratic differential zu we face the dilemma of either a formalization of each statement, 
with a precise formulation of the behavior of every trajectory in every situation, or 
a much more visual description, with rigorous proofs but illustrated by a number 
of figures. We opted for the second choice^ 

Next, we agree on some convention about trajectories. Let p(A g 7£ ■ be a zero of 
vj. We denote by 71 (p^)> 72 (p^)> • ■ ■, the trajectories of w emanating fr om p on 
7 Zj, in such a way that their canonical projections 77(7 n (p^))> see Section [T2I are 
enumerated in an anti-clockwise direction starting from the positive OX semiaxi^j 
Notice that when j s a branch point of 1Z, so that p ^ = p^ for some j ^ k, 
trajectories 7 n (p^) and 7 n (p^) are different because they emerge from p W = p^O 
on different sheets of 7 Z. Otherwise, when p ^ 1 belongs to a single sheet IZj, we 
occasionally drop the superindex (J) when it cannot lead us into confusion. 

Given two points p, q £ R, the integral 



along a contour 7 connecting p and q is understood to be the integral of any analytic 
continuation of the meromorphic differential \J—vo along 7. This is well defined up 
to the branch of the square root, which will be clear in each context. 


4.5.2. Degenerate case t = 0. For t = 0, the algebraic equation (3.9) reduces to 
(£ + z 2 ) (4C 2 - 8z 4 - 4£z 2 + 12z~3x 3 2/3 ) = 0, 


2 We confess we might have been influenced by the famous quote of Vladimir Arnold |82| : 

It is almost impossible for me to read contemporary mathematicians who, in¬ 
stead of saying “Petya washed his hands,” write simply: “There is a t\ <0 
such that the image of ti under the natural mapping t\ i—>• Petya(ti) belongs 
to the set of dirty hands, and a t 2 , t± < £2 ^ 0, such that the image of 1 2 under 
the above-mentioned mapping belongs to the complement of the set defined in 
the preceding sentence.” 

^ This notation is not correctly defined only if the direction of a trajectory coincides with a 
branch cut, situation that will be explicitly avoided in what follows. 
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whose solutions, denoted in accordance with ( 3 . 10 ), are 
y/ 3 y/h(z) 


Cl 0) = 


+ 


6(z) = --Z , &(/) = 


z 2 — a /3 y/h(z) 


2 w ’ ’ w 2 
where the branch of the square root is chosen to be positive for large real values, 

h(z) = 3 z 4 -4 z + 3 2 / 3 = 3 (z — a 2 ) (z - b 2 ) (z — b *) 2 , 
and the points a\, &i, a 2 , b 2l 6* are given explicitly by 


3 2/ 3 

ai = bi = —, 


b 2 = a 2 = + 


b * = 


1 

3V3' 


w 2 = 


The cut Ai is reduced to a single point a\, and the sheet 1Z - 2 is detached from 
the others. Since (3.9) is reducible, its Riemann surface is in fact the union of two 
Riemann surfaces, 

U 2 = C and 1Z = 1Z\ U IZ 3 . 

Here 1Z\ = IZ 3 = C \ A2, and A2 is a simple curve connecting the points a 2 ,b 2 , to 
be precisely specified later. The quadratic differential (3.11) degenerates into two 
quadratic differentials vj\ on C and w 2 on 7 Z, namely 
w\ = —3h(z)dz 2 on C, 

—1(3 z 2 — -\/3 ^h(z)) 2 dz 2 on7li, (4.22) 

—1(3 z 2 + y/3\/h(z)) 2 dz 2 on IZ 3 . 

We analyze the structure of their critical graphs next. 

Trajectories of zui, whose only critical points are as follows: 

• Simple zeros at z = a 2 ,b 2 ; 

• Double zero at 2 = 6*; 

• Pole of order 8 at z = 00. 

Under the change of variables z <—> z, the quadratic differential vj\ becomes 

— q(z)dz 2 , where 

q(z) = —( z 2 — 2 iz — 3)(z + i ) 2 

is, up to a multiplicative factor the same polynomial obtained in 27 eq. (2.1)]. 
Having in mind this identification, it was proven in 27 Theorem 2.1] that the 


trajectory 73(62) of w\ connects b 2 and a 2 : in the notation introduced above, 
73(62) = 71(02), see Figure [6j 

From [27 we also know that 73(62) intersects the real axis at a point a* < 6*, 
which can be calculated numerically: by 

a* « — 0 . 441782 . ( 4 . 23 ) 

The rest of the critical graph of w\ is as follows. Notice that due to the symmetry, 
we only need to describe the trajectories in the upper half plane. 

The trajectories 7/6*), 7/62), j = 1 , 2 , cannot be finite, see Principle P.3 above. 
Hence, they all diverge to 00 along the asymptotic directions ( 4 . 21 ), and according 
Theorem |B. 1 [ all directions are represented. 


There are 3 asymptotic directions 
for 4 trajectories in the upper half plane, so necessarily the divergence angle for 
71(6* *) is o[°°\ while the divergence angle for 72(62) is Since two consecutive 

trajectories emanating from a zero cannot diverge to 00 in the same direction (this 
would contradict Theorem B. 2 ), we conclude that both 72(6*) and 71(62) must 


diverge in the direction 9 ^, see Figure [bj 
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Trajectories of ?u 2 , whose only critical points are as follows: 

• Double zeros at z = z = bP and 2 = aP; 

• Double pole at z = oo* *- 1 ^ with real residue; 

• Pole of order 8 at oo^ 3 ). 

The double zeros bP = bp a 2 1 = lip are also branch points of 72, and the critical 
graph of W 2 is made of trajectories 'YjifiP), j = 1,2, emanating on 72.1, 

/o\ /o\ /o'N 

and of trajectories 7 j{a\ ), 7j(6 2 j = 1,2, along with 7 j(a\ '), j = 1, ..., 4, 
emanating on 7Z 3 , see Figure [7] 

The branch cut A 2 connecting the sheets 721 and 72-3, so far arbitrary, is chosen 
as A 2 = 73(62)) where 73(62) is the critical trajectory of w\ connecting 6 2 to a 2 , as 
described above. 

Lemma 4.6. With the branch cut A 2 specified above, the critical trajectories 
'Yjtyp) an d T j( a P) °f w 2 belong entirely to the sheet TZ\. 

Proof. Suppose that one of the trajectories Jj(bP) emanating from 6 2 intersects 
the cut A 2 for the fist time at a point x. Clearly, this point (actually, its canonical 
projection) must lie in the upper half plane: otherwise we readily get that r ij{bP) = 
7 i( a P) an d n0 intersection with A 2 occur. 

Integrating from 6 2 to x along 7 j(fiP) and using the definition of a trajectory 
we get 

fn px 

s 2 ds -— Re / ph(s) ± ds (4.24) 

where the ± sign in the last integrand depends on the side of the cut A 2 to which 
x belongs. However, A 2 projects onto the trajectory 73(62) of vo\, so the second 
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Figure 7. On the left it is represented the critical trajectories 
of n?2 on 771 and on the right its critical trajectories on 77-3, all 
of them in blue and just locally at the critical points. The black 
curve is the cut A2 connecting 771 and IZ 3 . 


integral in the right-hand side of (4.24) is purely imaginary. Hence, this equation 
reduces to 


Re ar = Re b';> = -. 


(4.25) 


It was proved in 27 that the trajectory A2 is contained in the domain bounded 
by the triangle with vertices 62, 3 -1 / 3 (v / 2 — 1) and a 2. In particular, the part 
of A2 on the upper half plane, and hence a;, is contained in the domain bounded 
by the triangle T determined by the vertices b 2 , Re b 2 and 3 _1 / 3 (v / 2 — 1). The 
function z 1 —> Re z 3 has a unique maximum on T at the point z = b 2 . Since Re z 3 is 
harmonic, it cannot attain a maximum on the domain bounded by T, hence (4.25) 
can only occur if x = b 2 , showing that 7j(&2^) \ {^} does not intersect the cut 

A 2 . “ □ 


Recall that 00 ^ is the center of a circle domain, which means that all trajectories 
of w 2 passing through sufficiently distant points on 77i are closed Jordan curves. In 
particular, no trajectory diverges to oo^ 1 ^, and every trajectory entirely contained 
in 72.1 has to be closed. Consequently, both trajectories 7,- (b^ 1 ) are closed as well, 
Jj(b^) = 7j(a2^)> an d 7i(^ 1) ) (respectively ^{b^)) intersects the real line, say 
at the point c» (respectively d*). We claim that 

d* < a* < c* < ai, (4-26) 


with a* defined in (4.23). 

Indeed, both c* and d* cannot lie on the same side of the cut A2 without running 
into contradiction with the general principle P.3 above. Hence, either d* < a* < c* 


or c* < a* < d*. The latter is impossible without 71 (&2 1 '*) an, 4 72 (b^) intersecting 
somewhere in the upper half plane, which again contradicts P.3. We conclude that 
d* ^ cq ^ c*. 
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Let us prove the inequality c* < ai. Using the definition of trajectory and (4.19) 
we get 


r(U 


0 = Re 




Q(s)ds = Re / (&(«) - &(s))ds = w 3 + / (£ 2 (s) - £ 3 (s))ds. 


' b 2 


But w 3 = w 3 (0) < 0, see Figure |4j so that 

[ (6(s) - &{s))ds > 0. 


Function £2 — £3 is continuous and non-vanishing on (a*,+ 00 ), and by (3.10) it is 
negative for large real parameters, so it is negative on the whole interval (a*, + 00 ). 
Since c* £ (a*, + 00 ), the equality above is only possible if c* < a±. This proves 

(p6). 




Figure 8. r = 0: blue lines represent the critical graph of 1372 on 
72 1 (left) and 72 3 . The black curve is the branch cut A 2 connecting 
72 1 and 72 3 . 


The discussion of the structure of the trajectories 71 (a^), 72 (a^), hib^), 
72(^2 ) is identical to the analysis of the trajectories 77j(& 2 ), j = 1 , 2, on 72 2 
for 737 i above, so we omit the details. 

The global structure of the critical graph on t 37 2 on both sheets is presented in 
Figure [HJ The basic conclusion is that with the branch cut A 2 specified above, the 
critical graph splits into two sets: a closed Jordan curve on 72. 1 , containing and 
? 4 ' 1 , and 4 analytic arcs on 72 3 , starting and ending at oo^ 3 \ each passing through 
one of the branch points ab ^ and a[ 3 K 

4.5.3. Trajectories for 0 < r < 72 • A combination of the general principles P.2 and 
P.6 assures us that the finite critical trajectories for r = 0 remain finite for small 
perturbations of t, and that the behavior of the trajectories described for r = 0 is 
preserved for r small. 

Let e > 0 be sufficiently small. The general principle P.4 above tells us that if 
we consider the domains (respectively and fi^), swept by trajectories 
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Figure 9. Left: critical graph of vo for r = 0, with the domains 
j = 1,2,3, in gray. Right: local behavior or the critical 
trajectories for vo and r = 6 > 0 , passing through a^ \ j = 1 , 2 ,3, 
with the same domains superimposed. 


of (4.221 passing through points in the e-neighborhood of (respectively a{ ^ 
and dj 3 ^), then there exists a S > 0 such that the critical trajectories for vo and 
0 < r < <5, passing through a^ \ j = 1,2,3, belong to fl E = U U fl E 3) . 
These domains are depicted schematically on Figure [9] 

We should keep in mind that the branch cut A 2 now is completely specified by 
which is consistent as long as the critical trajectory 73(62 j joining 
exists and remains on 1Z 2 for the full range of the parameter r under 


Definition 
ai 2 ' 1 and b. 


4.5 

(27 


consideration. The forthcoming analysis shows that this is the case for r < r c . 
According to the general principle P.4, there exists a S > 0 small enough such 


that for 0 < r < 5 , both a[^. b^ 1 ' 1 € QW, j = 1,2,3, and the critical trajectories 
emerging from af\b^ stay in 


,(i) 
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For instance, since for r = 0 the critical trajectories 7 j = 1,..., 4, 
define three half-plane domains (bounded by 7 i(a^) 1172 ( 12 ^), 73 (a 3 3 ^) U 74 (a^) 
and 74 ( 0 ^) U 74 (a^), see Figure |ij right), they must be persistent under small 


/ o\ 

perturbation of r, and either a) or b'f 1 , or both, must belong to their boundaries. 

(3) 

Taking into account the structure of it is straightforward to conclude that for 
0 < r < S, the trajectories of w through a( Jj> and b < (‘ > are as shown in Figure 


l(3) 


We now examine the trajectories from and b\ x> . The following result comes 
in very handy: 




10 


Proposition 4.7. Let 


h(x,y) := 


I Re(£i + (s)-£ 3 (s))ds= f Re(£ 2 +( s ) ~{ 3 (s))ds, if x,y 

J X J X 


e Ai, 


f Re(£i(s) — ^ 2 +( s ))ds = f Re(£i(s )-£ 3 +(s))ds, if x,y € A 3 , 

*• J X J X 

(4.27) 

where we integrate along each interval. 

If 0 < t < 1/12, then there exists no pair of values 1 / 1 /, x,y £ Ai, such that 
h(x, y) = 0 . 

If 1/12 < r < 1/4, and there exists a pair of values x 7 ^ y, x, y £ Ai (resp., 
x,y £ A 3 / such that h(x,y) = 0, then there exists no such pair of values on A 3 
(resp., on A\). 

Furthermore, if x 1 < 2 / 1 , X 2 < 2/2 are two such pairs, ( 27 , 2 / 1 ) D ( 27 ) 2 / 2 ) 7 ^ 0- 


Notice that h is well defined on A 1 and A 3 due to the symmetry relations (4.9) 
and (4.12|, and that A 3 = 0 for 1/12 < r < r c . 

Proof. Assume that there does exist a pair of values x < y, 2 :,y £ Ai, such that 
h(x,y) = 0 (same analysis is valid for x,y £ A 3 ). By the mean value theorem, 
there exists a u £ (x, y) such that 


hence, 


But 


Re£i+(M) = R e( 2 +(«) = &(«), 

0 = Re£i+(u) +Re£ 2 +(«) +£ 3 (u) = 3£ 3 (u). 


D{u) = + {u)i 2 +{u)^{u) = 0, 

and the assertion follows from Proposition [42] keeping in mind that D has exactly 
one zero on A 3 U A 3 for 1/12 < r < 1/4. □ 

One of the consequences of Proposition [47] is that for 0 < r < <5, the trajectories 
emanating from a 3 1 ' 1 and b^ cannot cut A 3 and must stay on the sheet IZ\. Thus, 
using again the general principle P.4 we conclude that the trajectories 7 / c (< 4 1) ), 


10 


7 fc(&i ), k = 1,2, are closed and encircle the cut A 2 , see Figure 

Similar considerations can be applied to get the behavior for the trajectories 
emanating from a^\ b^\ and the final result for r small is seen in Figure [lo| We 
skip the details. 

The outcome of our analysis is that the critical graph of zu has the structure 
showed in Figure [lOl at least for 0 < r < S. Our next goal is to prove that this is 
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Figure 10. Critical graph of vo for 0 < r < 1/12, with the strip 
and ring domains labeled by S/’s. Notice that some of these do¬ 
mains intersect more than one sheet. 


actually valid for r £ (0,1/12). The continuity principle P.4 yields that this is the 
case as long as 

(i) No collision of the critical points occur: this is true indeed for r £ (0,1/12), 
see Section 14741 

(ii) No new domains emerge, which amounts to say that the finite trajectories for 
0 < t < S remain critical for 0 < r < 1/12: this is assured by a combination 

of P.2 and P.6. 

(iii) No connected components of the complement of the critical graph “disappear”. 
More precisely, it means that no width of any strip or ring domain becomes 
zero. These domains for 0 < r < S are identified on Figure [lOj there is one 
ring domain S± and 4 strip domains, S 2 , ■ ■ ■, S 5 . There widths cr(Sj) (see the 
definition (B.6) in Section [B.3| are given by: 

• cr(Sj) = \ojj |, j = 1,2,3, as defined in (4.19). They do not vanish for 
r £ (0,1/12), see Figure [IJ 

• cr(S , 4 ) = |/i(ai,6i)|, with h defined in (4.27), which does not vanish for 
t € (0,1/12), see Proposition |4.7| 


a ( S 5) = f b ‘(£ l(s) - &s(s))ds 

E6 


see 


which does not vanish for r G ( 0 , 1 / 12 ), 


We conclude that the critical graph of m, depicted in Figure |lT)| is valid for the 

‘ / o \ 

whole range 0 < r < 1/12. In particular, the critical trajectory 71(02 connects 
(2) (2) 

the points Oj and b\ , which proves Theorem 


4.4 


for 0 < r < 1 / 12 . 


4.5.4. Trajectories for ^ < r < t\. When r = 1/12, the double point 6 * coincides 
with bi, and the strip domain S 5 disappears (<j(S 5 ) \ 0 as r /• 1 / 12 ), see Figure |TTj 
left. Clearly, this transition has no impact on the structure of trajectories on the 
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Figure 11. Left: critical graph of ro for r = 1/12, with the do- 
mains 12/, j = 1,2, in gray. Right: local behavior or the critical 
trajectories for w and r = 1/12 + 5 > 0, passing through b\ 3 \ 
j = 1, 2, and b^\ with the same domains superimposed. 


third sheet. Moreover, again a combination of P.2 and P.6 assures the finite 
trajectories for r = 1/12 remain finite for 1/12 < r < 1/12 + 5. 

Let e > 0 be sufficiently small. Similarly to what has been done in the previous 

interval, the general principle P.4 tells us that if we consider the domains $2^ and 
( 2 ) 

12 / , swept by trajectories of m passing through points in the e-neighborhood of 
and b[^ = b^\ then there exists a 5 > 0 such that the critical trajectories 
for 1/12 < t < 1/12 + 5, passing through b± \ j = 1,2, and b^\ belong to 
12 e = U 12^. These domains are also depicted schematically on Figure fllj 
left. 

For 1/12 < r < 1/12+5 we consider the first sheet and the trajectories emanating 
from b^ and b^; thanks to principle P.2, we concentrate on the upper half plane 
C + (or to be precise, on its pre-image by 7 r on TZi), namely 71 72 (5*^), 
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7 1 (6^), see Figure [TTJ right. These trajectories must stay in so they have to 
intersect 7 r _ 1 (K) on IZi. Let us denote the points of intersection of 7 i( 67 ), 72( 67 ), 
71 (& 7 ) by x^\ respectively. Using the general principles P.2 and P.3 

we must immediately discard the following possibilities: (i) Xj > b± for some j, 
(ii) < a 7 and X2 < a 7 (recall that a* = A 2 flK). Since trajectories cannot 
intersect, it holds x\ < x 2 < 2:3 and we conclude that necessarily x 2 £ (oi,&i), 
and consequently, X3 € (ai,6i) as well. In particular, 6(23,61) = 0 , in the notation 
( | 4 ~ 27 ). 

Since X\ and 2 2 belong to trajectories with a common point & 7 > the assumption 
x\ £ (ai,6i) yields that 6,(21,22) = 0, and since (xi,x 2 ) D (2:3,61) = 0, this 
contradicts Proposition | 4 . 7 | 

From the considerations above, it follows that 71 (& 7 ) is closed, stays on IZi, and 
intersects 7 t _ 1 (R) to the left of A 2 , and the trajectories 72( 67 ), 71 (& 7 ) intersect 
the cut Ai and move to the second sheet 1Z 2 - We keep denoting these points of 
intersection by x 2 ,X3 as before. Clearly, 2:3 > 2 2 , and Proposition | 4 . 7 | implies 
that these are the only points of intersection of these trajectories with the interval 
(«i, 61 ). 


/o\ 

Let us turn to the second sheet, 1Z 2 and consider 71(6) ): from the structure 
of fl7 it is clear that it either diverges to cxd 2 \ or intersects the branch cut 
between ai and 61 and moves to 1Z\. If we assume the latter, fl 7 shows that 
either it will return to TZ 2 at a different point in (ai,&i) (in contradiction with 
Proposition 4 . 7 ), or it bounds a simply connected domain, which contradicts P.3. 
Hence, 7i (& 7 ) must diverge to oo^ 2 \ and thus &7 lies on the boundary of the half 


(n\ (2) 

plane domain bounded for r = 1/12 by 71(61 ) and 74 ( 6 ) ). In particular, by P.5, 
for 1/12 < r < 1/12 + 6 , 7 i(& 7 ) diverges to oc/ 2 ) in the same asymptotic direction 
given by the angle from ( |4.21 1 . 

Recall that we concluded that the trajectory 71 (&7) enters 1Z 2 through the cut 
(ai, 6 i) a t a p oint X 3 . The only possibility left for it is to go to oo^ 2 b Applying 
Theorem B.2 to the tu-polygon bounded by the trajectories 7 i( 6 ^) and 72 ( 6 ^); 
we get that 71 (& 7 ) g° es to oo*- 2 ^ in the asymptotic direction given by the angle 
an b consequently 71 ( 67 ) extends to cx/ 2 ) with angle d 5 as well. 

The outcome of our analysis is that for 1/12 < r < 1/12 + 6 the critical graph of 
vo has the structure showed in Figure [12] We prove that this is actually valid for 
t £ (1/12, n). Again, the continuity principle P.4 yields that this is the case as long 
as (i) no collision of the critical points occur (this is true indeed for r £ (1/12,1/4), 
see Section 4.4), (ii) finite critical trajectories for 1/12 < r < 1/12 + 6 remain finite 
for 1/12 < r < ti (assured by a combination of P.2 and P.6); and (iii) no width of 
any strip and ring domains become zero. These domains for 1/12 < r < 1/12+6 are 
identified on Figure fl 2 | there is one ring domain S± and 4 strip domains, S 2 , ■ ■ ■, S 5 . 
There widths v(Sj) are given by: 


cr(Sj ) = | ujj\, j = 1,2,3, as defined in (4.19). They do not vanish for 
1/12 < r < Ti, see Figure [4] although cr(S 1 ) does vanish for t = T\. 
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Figure 12. Critical graph of w for 1/12 < r < ri, with the 
strip and ring domains labeled by Sj s. Notice that some of these 
domains intersect more than one sheet. 


Re 


is given by the absolute value of 

[ ( 66 ) - 66 )) = Re f ( 66 ) -&(s))ds + Re f ( 66 ) - 66 )) ds 
/ a\ J a\ J b<2 

— —UJ± + CJ4, 

and wi 0 W 4 for r > 1/12, see Figure [4] in Section 4.3 


(4.28) 


• <t(<Ss) = / hi *(60) —6 {s))ds , which does not vanish for r > 1 / 12 , see 

We conclude that the critical graph of w, depicted in Figu re [I2l is valid for the 
whole range 1/12 < r < ti. In particular this yields Theorem |4.4| in the mentioned 
range of r. 


4.5.5. Trajectories for t\ < t < r c . At the value r = t± the critical trajectory 
7162 ^) on the first sheet hits the branch point < 2 ^, so that the ring domain S/ 
disappears (cr(5'i) \ 0 as r / r/, see Figure [12} From the analysis of the behavior 
of the rest of the widths cr(Sj) and the other finite critical trajectories, it follows 
that this fact does not affect the rest of the strip domains - note that there are no 
other ring domains. In particular, the trajectories emerging from ai, and b 2 on 
the sheets 1Z 2 , IZ 3 do not display any phase transition. 

The critical graph for r = t\ is depicted in Figure |13| In accordance with 
the methodology we have followed so far, we fix an e > 0 sufficiently small and 
consider the domain Q e swept by trajectories of vj passing through points in the 
^-neighborhood of 1 and a ^ . Notice that no longer lives on the single 

sheet, and its boundary now also contains critical trajectories, namely 7 j(a) ), 
7 j( 4 3) ), 7 j(^ 2 3 ^)? J = 1)2. This is so because, as already mentioned, there is no 
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Figure 13. Left: critical graph of vo for r = ti, with the domain 
f l e in gray. Right: local behavior or the critical trajectories for w 
and r = n + S > 0 , passing through b^ and a^, with the same 
domain superimposed. 


transition for these trajectories. Observe also that when r /• t\ the trajectory 
72 = 72 (<4'^) does not collide with any critical point other than its endpoints, 

hence its topology is unchanged under small perturbations of r around t\. In fact, 
there exists a S > 0 such that the critical trajectories for T\ < r < 17 + 5, pas sing 
through b^\ belong to f l e . This domain is also depicted schematically on Figure 13 
left. 

We now consider the possible behavior of 71 (b^) for t\ < r < 77 + <5, having in 
mind that it cannot leave the shaded region which shows that either 71 (b^) 
moves immediately to 7 Z 3 through A 2 , or it intersects the (preimage of) the real 
line near a 1 . 

In the first case, 7 i(& 2 1 ^) extends to oo^ 3 ^ with angle Thus 71 
is the boundary of a zn-polygon for which k = 1 and A = 0, contradicting (B.7|. 
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Figure 14. Some trajectories of zu on the first two sheets, for 

Tl < t < T c . 


Assume otherwise, so that 71 ( 62 ^) intersects the (preimage of the) real line at 
a point irW, x close to a 3 . If x < a±, we integrate from bo to x over 71 (b 2 ^) and 
then from x to ci\ over the real line to get 

nX na\ na\ 

wi = Re / ( 6 (s)- 6 (s))ds+ / ( 6 (s)- 6 (s))ds = / (&(s)-£, 3 {s))ds. (4.29) 

J 62 J X */ X 

In the range of r considered, > 0, while the last integral is < 0, see ( |4.5| ), which 
leads us into a contradiction. Hence, x £ (ai,&i), so that 71 (f*^) moves to the 
second sheet. Recall that if £ 3 ^ is the point o f int ersection of 72 (f^ 1 ') with Ai, we 

7i (^ 2 ^) cannot return to 7Zi; in 


4.7 


already have h(x 3 ,bi) = 0, so by Proposition 
consequence, it stays on the second sheet and diverges to oo^ in the asymptotic 


14 


direction corresponding to the angle see Figure 

On the other hand, trajectories 7 Aa^), j = 1,2, are “trapped” between 73 (b^) 
( 2 ) 

and the branch cut Ao , see Figure 14 left. Thus, they cannot stay on 77 1 without 


violating the general principle P.3, so they move to the third sheet and diverge 
to 00 in the asymptotic directions corresponding to the angles fig 00 '* and 8 ^°\ 
respectively. 

The outcome of our analysis is that for n < r < t\ + S the critical graph of w 
has the structure showed in Figure [l5j Same analysis as in the previous section 
shows that this is actually valid for the whole range 77 < r < r c , in particular 
implying Theorem |4.4| for this range of r. 

4.5.6. Trajectories for t c < r < T 2 . We come to the topologically most important 


phase transition. According to Definition 4.5 we use the lift of the trajectory 


joining a' 2 and b y 2 ’ on 1Z 2 as the branch cut connecting the sheets IZi and 1 Z 3 . 

When r = r c , the strip domain S 2 disappears (cr^) \ 0 as r t c ), and this 

( 2 ) 

trajectory finally reaches the branch point a) , see Figure 


15 


Note that this 

transition corresponds to (|1.34|). We fix an e > 0 sufficiently small and consider the 


domain fL swept by trajectories of w passing through points in the e-neighborhood 
of a\ , a 2 , and b 2 . The critical graph for r = t c along with fl e is displayed in 
Figure [l6| 

For t c < t < t c + S we examine first the second sheet and the trajectories 
emanating from a ( j 2 '*, a 2 ' 1 . and b^\ For a sufficiently small S > 0, these trajectories 
must stay in see Figure [l7| 
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Figure 15. Global structure of trajectories for n < t < t c . 







Figure 16. Critical graph of vo for r = r c , with the domain f l £ 
in gray, consisting of the trajectories for w passing through the 
e-neighborhood of a^\ and ■ 


(O') 

Analyzing 71 (a) ) we have to discard the following possibilities: 

• 71 ( 0 ^) cannot intersect (the preimage by tt of) the real line to the left of 
ai without violating the general principle P.3. 
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Figure 17. The domain fl e on 72-2 for t c < t < t c + 5. 


it cannot intersect (the preimage by 7 r of) the real line to on the cut (ai, 6 i) 
either: otherwise the equation h(x,a\) = 0 has a solution in x £ (ai,&i), 
along with the identity h(x 3 , bi) = 0 , where ig 1 "* is the point of intersection 
of 72 ( 6 ^) with Ai, and since (a 3 ,x) fl ( 0 : 3 , 61 ) = 0 , this would contradict 
Proposition |4.7| 

7 i(a^ 21 ) 73(62 ), because 1 x 2 ^ 0 for r > r c , see Figure [ 4 J 

7 i(a^ 2 ') cannot diverge to 00 ^ in the asymptotic direction d ^ 00 ' 1 ■ Indeed, 
otherwise either 73 ^J^ 2 "*) also diverges to 00 ^ in the same direction, or it 
intersects the real axis to the left of a\ . In the former case, we get a txj- 


polygon for which n = 1 and A = 0, contradicting (B.7), and in the latter 


one we proceed as in (4.29) (with uq replaced by uj 2 ) to get a contradiction. 

The only possibility left for 71 (a^) is to diverge to 00 ^ in the asymptotic 
direction . Since it was already observed that 73 ( 6 ^) cannot diverge to ocA 2 \ 
it must intersect 7 r( _ 1 l(R) to the right of ai. 

The outcome of our analysis on the second sheet, as well as the region f 1 e on the 
remaining sheets, is displayed in Figure |18| The cut A 2 is chosen in such a way 
that its projection on 7l 2 coincides with \^y 3 {b^) U 71 ( 02 ^)^ (11Z 2 . 

What is left is to describe the behavior of 73(63 ; ) on the rest of the sheets. 
We already saw that this trajectory has to move to IZi through the cut A-]. It 


cannot intersect A i ag ain (see Proposition 4.7), and it must stay in the region f l E 
displayed in Figure 18 Hence, the only possibility is that 73 ( 6 , ) intersects the cut 
A 2 , moves to the sheet 7^3 and diverges to oo^ 3 ) in the asymptotic direction 6 4 °°^. 

The critical graph of w for r c < r < t c + 8 has the structure showed in Figure [l9l 
and we now prove that this is actually valid for r £ (t c ,t 2 ). Again, the continuity 
principle P.4 yields that this is the case as long as (i) no collision of the critical 
points occur (this is true indeed for r £ (t c ,T 2 ), see Section |4~4| ); (ii) finite trajec¬ 
tories for t c < t < t c + 5 remain finite for r £ (t c ,t 2 ) (assured by a combination of 
P.2 and P.6) and (iii) no width of any strip domains become zero. These domains 
for t c < t < t c T S are identified on Figure |19| observe that there are now 7 strip 
domains Sj, and no ring domains. Let us compute their widths cr(Sj). 

Lemma 4.8. For t c < r < r 2 , 


tr(S'i) = <t(S 2 ) = <j(S 3 ) = |wi(t)|, 


with uj(t) defined in (4. 19). 
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Figure 18. Part of the critical graph of vo for t c < t < t c + 5, 

with the domain fin gray, consisting of the trajectories for w 

(2) (2) (2) 

passing through the e -neighborhood of a\ , a 2 and 61, • 



Hi 
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Figure 19. Critical graph of w for r c < r < r 2 , with the strip 
domains labeled by Sj’s. 


Proof. The fact that cr(Si) = |wi| is the straightforward consequence of the defini¬ 
tion of uj\ . Also <j(S-2) = cr^Ss) by the symmetry under conjugation. So, it remains 
to compute the width cr^), given by the absolute value of the real part of the 
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U 2 ) 





Q(s)ds ■ 




Q(s)ds 


pa i ra<2 

/ (Ci {s)~&(s))ds+ (Ci(s )-6 (s))ds. 

J 62 J a\ 


Symmetry under conjugation tells us 


pCL2 pb2 

Re / (Ci(s) -f 2 (s))ds = Re / (£i(s) - &(s))ds, 

«/ Cli «/ Gti 

and using it in the previous identity, we get 


cr{S 2 ) = 


Re 


pa 1 pa\ 

/ (Ci(s) - &(s))ds - Re / (Ci(s) - C 2 (s))ds 

J bo J bo 


I b 2 
r a 1 


Re/ (&(s) - 6 (s))ds 
Jb 2 


b 2 

= |wi|. 


□ 


Lemma 4.8 implies that cr(Sj) ^ 0, j = 1,2,3, for r c < r < r 2 , although 
cr(S'j) \ 0 as t t 2i j = 1, 2, 3, see Figure |4j 
Regarding the rest of the strip domains, 

• cr(S±) = |oJ 4 1 7 ^ 0 for r c < t < r 2 , see again Figure [4] 

• 0 / 55 ) = / h h *(6(s) ~ C 3 {s))ds , which does not vanish for r > 1/12, see 

flirt . 

• From the symmetry by complex conjugation, a(Se) = <r(S 7 ), and the struc¬ 
ture of 56 on 7*1 2 shows that ct( 56) = |w 2 |, so that cr(Se) = cr(Sr) 7 ^ 0 for 
t c <t < t 2 . 

Consequently, the critical graph displayed in Figure [19] is valid for r c < r < t 2 . 


In particular it implies Theorem 4.4 for r c £ (t c ,t 2 ). 


4.5.7. Trajectories for t 2 < t < 1/4. Lemma 4.8 shows that at r = r 2 , the strip 
domains 5i, 5 2 and 53 displayed in Figure m disappear simultaneously , which 
happens because at that moment the branch point hits the critical trajectory 
7 i(Q' 9 1 ') = 72 (fr? 1 *) and the point a* of intersection of A 2 with the real line collides 
with the critical trajectories 72 (a^) and 71 (fr^). The resulting critical graph is 


shown in Figure 20 


As before, we fix an e > 0 sufficiently small and consider the domain swept 
by trajectories of w passing through points in the e- neighborhood of ckf*, b^ 
and a* \ Observe that for a small perturbation of r = r 2 , 73 ( b 0 7 Z 2 does not 
coalesce with critical points other than the starting point b 2 , hence this arc of 
critical trajectory displays the same structure as for r = r 2 . In particular, the cut 
A 2 is well defined as the projection of (73 (b\ ) U 71(03 ')) 0 7*1 2 on the other sheets. 
Thus, there exists a <5 > 0 such that the critical trajectories for r 2 < t < r 2 + <5, 
passing through the above mentioned critical points, belong to Cl e . This domain is 
also depicted schematically on Figure [20j left. 
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Figure 20. Left: critical graph of vo for r = T 2 , with the domain 
in gray. Right: local behavior or the critical trajectories for w 
and t = T2 + 6 > 0 , passing through a^\ and a*, with 
the same domain superimposed. 


Lemma 4.9. Let t be a point on the part of the curve joining b^ with ai 1 '*, and 
7 a Jordan curve on 1Z\ connecting the boundary values t± on and containing 
the only branch point blp inside. Then 
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Proof. We can deform 7 to the cut A 2 . If we denote by A 2 (f) the arc of A 2 from 
b 2 to t, we get 


Re [ Q(s)ds = Re [ (£ 2 +(s) - £ 3 + (s))ds - Re [ (£ 2 _(s) - £ 3 -(s))ds 

J 7 J A2(t) J A 2 (£) 

= Re [ ( 6 s- 0 ) - 6+(s))ds = Re [ (Ci +(s) - ^ 3+ (s))ds 

J A 2 (t) J A2 (t) 

i / Q + ds = 0 , 


= Re 


1^—1 /2) 

where we have used the jump condition ( 4 . 11 ) and the fact that AJ, ; is an arc of 
trajectory of tu. □ 

Now we describe the critical trajectories for r 2 < r < t 2 + <5, starting with 
72(a) 1 )). Notice first that in this case uj± 7^ 0 , so that y 2 (a) 1 )) cannot contain 
&) 1 ). Furthermore, if 72(a) 1 )) intersects the cut A 2 (and thus diverges to oo( 2 ( 
in the asymptotic direction then 71(a) 1 )) U 72(a) 1 )) U 7i(&) 3 )) determines a 

cu-polygon for which k = 1 and A = 0 , in a contradiction with (B. 7 ). 

Keeping in mind that 72(a) 1 )) must belong to f l e we conclude that it has to 
intersect the real axis in one of the sheets. Let us denote by x the point of the first 
intersection of this trajectory with 7r _1 (M). Again, we discard some cases: 

• x cannot be on 77 .1 to the left of ai: this yields (by the general principle 


P. 2 ) that 7 2 (a) 1) ) = 73(a) 1,1 ), violating P. 3 . 
x cannot be on 77 .1 to the right of a*. 

Indeed, if x > a*, we form a curve 7 given by the union of three pieces: 




7! is the arc of 71(a) 1 )) from a* 1J to a) l; , 7 2 is the arc of 72 (a) i; ) from 
to x, and 73 is the interval from x to a) 1 ). This curve satisfies the 


,(i) 


(i) 


AW 


,(i) 


assumptions of Lemma 4.9 so that 


R-e [ (6+(s) - &+{s))ds = 0. 

J -V 


Since 


Re 


f (6+(«) - &+{s))ds = 0 
J 71-U72 


by the definition of trajectories, we conclude that 

Re f (£ 2+ (s) - £,3+{s))ds = 0 , 

J a * 


that is, h.(a*,x) = 0 , in the notation ( 4 . 27 ). But if a;) 1 ) is the point of 
intersection of 72(6^) with A l5 we also have h(x 3 ,bi) = 0 , and since 
(a*, a;) (7 (x 3 ,bi) = 0 , this would contradict Proposition 4.7 


x cannot be on 1 Z 2 to the left of a* (in other words, x cannot belong to A 3 , 
which could occur if the trajectory 72(a) 1 )) had slipped to the second sheet 
through the cut A 2 before hitting the real line): this yields 


Re 


(£1 +(s) - f 3 +(s))ds = 0, 
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S 5 


'■--..S 4 \s 2 \ Se 

Si / S 6 

.t. 

S 5 





Figure 21. Critical graph of vo for T 2 < r < 1/4, with the strip 
domains labeled by Sj’s. 


„(i) 


that is, h(x,a *) = 0, in the notation (4.27), with x,a* £ A 3 . But if £3 ' is 
the point of intersection of 72(^1 ) with Ai, we also have h(x 3 , 61 ) = 0 , in 
contradiction with Proposition |4.7[ 


The only possibility left is that 72 (<4^) hits the real line precisely at the point 
a*. This means also that 72 ( 0 ^) = 71 ( 0 ^), as shown in Figure 21 

Regarding the trajectories 71 (& 2 ^) an d 72 ( b 2 ^), they cannot satisfy 71 (fr^) = 


72 (^ 2^)1 an d due to principle P.3 they must belong to f l £ . Hence, they have to 
behave as shown in Figure [2lj We skip the details. 

The critical graph of w for r 2 < r < r 2 +<5 has the structure showed in Figure [2l| 
and we prove that this is actually valid for r £ (t 2 , 1/4). The continuity principle 
P.4 yields that this is the case as long as (i) no collision of the critical points 
occur (this is true indeed for r £ (72, 1/4), see Section [OJ, (ii) the finite critical 
trajectories for T 2 < r < 72 + 6 remain finite in the range r £ ( 72 ,1/4) (certainly 
true by a combination of P.2 and P.6) and (iii) no width of any strip domains 
become zero. These domains for 77 < r < 72 + S are identified on Figure |21| 
observe that there are now 6 strip domains Sj, and no ring domains. Their widths 
cr(Sj) are: 


cr(S , 1 ) = a (S 2 ) = |wi| ^ 0 for r £ (r 2 ,1/4). 

0 /S 3 ) = a(S 4 ) = |w 3 j ^ 0 for r £ (r 2 ,1/4). 

a/S/) = / fc h /(^ 2 (s) ^£3 {s))ds , which does not vanish for r > 1 / 12 , see 

fOT). 


As in (4.28), we get that a(S 6 ) = \u 4 — Wi|, and 07 7 ^ uj 4 for r > 1 / 12 , see 
Figure |4j in Section |4~3j 
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We conclude that the structure of the critical graph of w, depicted in Figure [21] 
is actually valid for the whole range t £ (r 2 ,1/4). This finishes the proof of Theo¬ 
rem mu 

Remark 4.10. The attentive reader might notice that in Figure [l9j for instance, 
the critical trajectories 7 i(c 4 2 ' > ) and 72 (a] 1 ' 1 ) intersect the cut A 2 in IZi on pairs 


of opposite points t±. This phenomenon, which also occurs in Figure 21 is easily 


explained by Lemma 4.9 


Remark 4.11. It follows from the results of Sections 4.5.3 -4.5.7 that the trajectories 
(2) (o) 

7i(& 2 0 and 72(^2 ) determine a half plane domain H. In particular, this implies 
that there is an orthogonal critical trajectory 71 emerging from b 2 ; which is entirely 
contained in H and extends to oo *- 2 - 1 along the angle 27r/3. Similarly, there is an 
orthogonal critical trajectory 72 (which is the complex conjugate of 71 ) emerging 
from a 2 ~ ] and extending to 00 ^ along the angle —2tt/3. Then the projected contour 

r = 77(7! u 72) u a 2 


satisfies the conditions stated in (3.13) (3.15). 


Appendix A. Proof of Theorem 11.31 in the general case 

presented in Section [ 2 ] is valid when the set S a C S Q 


The proof of Theorem 


1.3 


of double poles of the coefficient D in (1.13) is empty (see Proposition 2.8). Our 
goal now is to show that for a vector of measures p £ A4 a whose components are 
supported on a finite union of analytic arcs, and such the associated functions 


in (1.17) satisfy (1.13) for a polynomial R and a rational function D , the equality 

T> h (p)=0 (A.l) 

is valid for every function h £ C 2 ( C), without any further restriction on the poles 
of D. Recall that this fact was established so far for the Cauchy kernels h z defined 
in (12.31), and when supp hf] S a = 0, see Proposition 


As a first step, we extend (A.l) to polynomials: 


2.1 


Lemma A.l. Under the assumptions of Theorem 1.3, P q {p) = 0 for every alge¬ 
braic polynomial q. 


Proof. Fix p > 0 for which 

supp pi U supp P 2 U supp p 3 C V p := {.t £ 
For x £ V p and z sufficiently large, we can expand 

r.i 


\x\ < p). 


h z(x) = ~Y1 


Z3 + 1 ’ 


1=0 

which implies the identity 

p n -i{x, z) - z n+1 h z {x) = x n + w z (x), 

where the function 


(A.2) 


•Xx)= Y, 


3—n+l 


Z3~ 
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converges to 0 as z —>• oo uniformly for x £ V p , and p n -i(x,z) is a polynomial of 
degree n — 1 in i, given explicitly by 

n— 1 

Pn-l(x,z) = - T; X 3 Z n ~ 3 . 

3=0 

Set p_i =0. If 


N 


a n x 


(A.3) 


q(x) = ' 

n—0 

is a polynomial of degree N, then ( fA~ 2 | ) shows that 

g(x) = ~zq(z)h z (x) + Q z {x) + W z (a;), 

where 

TV AT 

Qz{x) = y: anPn-i(x, z), W z {x) = - ^ a n w z (x). 

n =0 n—0 

Note that Q is a polynomial of degree N — 1 in x. 

Moreover, W z and W' z both converge to 0 as z —> oo uniformly for x £ V p . Hence, 
the convergences 

W.(z) - W,(y) ^ 0 V(x)Wz(x) ^ 0 

x-y 

hold uniformly for x,y £ V pi and from the definition of T>h in (2.1) we get 

V Wz (ft) ~> 0 as z —> oo. 

The quantity Vh(jl) is linear in h , so ( |A.2 ) implies 

Vqijl) = -zq{z)V hz (p) + V Qz (p) + V Wz (p) = V Qz (p) + V Wz (p), 

where for the last equality we used Corollary |2.7[ and hence from the previous limit 
we get 

V q (p) = lim V Qz (p). 

z—>oo 

The result now follows easily by induction on the degree N of q. If N = 0, then 
Q z is identically zero, and the equality above implies that 'Dq(p) = 0. Assuming 
now that T>h{p) = 0 for every polynomial h of degree at most N — 1, we get that 
T>Q z (jl) = 0 because Q z {x) is a polynomial of degree at most N — 1 in x, and the 
equality above implies 'Dq(p) = 0 , concluding the proof. □ 


If p £ S a , then as it is discussed in Remark 


1.10 


the point p belongs to exactly 


two of the supports of pi, P 2 and p 3 , and locally the union of these sets is an 
analytic arc. That is, there exists an open disk U p centered at p such that 

7p = Up fl (supp pi U supp p 2 U supp p 3 ) (A.4) 

is an analytic arc passing through p. 

Additionally, given any point 2 £ (supp pi U supp supp p 3 )\S a , there exists 
a small disk B z centered z, disjoint from S a and such that the set 

B z n (supp pi U supp p 2 u supp p 3 ) 


is a finite union of analytic arcs, which can only intersect at the common point 2 . 
In case z £ S a , this intersection reduces to a single analytic arc. The collection 

{B z } U {U p } 
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constructed above is an open cover of the compact set supp fi i U supp /i 2 U supp /u. 3 , 
from which we extract a finite subcover 


{Bj} j=1 U {U p } p& § p , 

where for j = 1,..., m we have Uj = U z for some z = Zj £ (supp /ii U supp /i 2 U 
SUPP/Z 3 ) \ S a . Set 

m 

U= (J U p . 

J =1 p&s a 

It follows from their construction that these sets satisfy 

B n S a = 0, supp Hi U supp fi 2 U supp C B U U. (A.5) 

Consider a smooth partition of unity {ipk} of supp/xi U supp/r 2 U supp ^3 sub¬ 
ordinated to the open cover B U U. That is, each function ipk is real, belongs to 
C'oo(C) and additionally satisfies the following properties. 

• 0 < ipk(z) < 1) for every z £ C. 

• For every k, supp 1 / 7 , C U or supp^fe C B. 

• Every z £ supppi U supp ^2 U supp/ 1,3 belongs to the support of a finite 
number of functions in the collection {ipk}. 

• Sfc i’k(z) = 1 , for every 2 £ supp/ii Usupp/r 2 Usupp/x 3 . 

Since supp /Ti Usupp p 2 Usupp H 3 is compact, the collection {ipk} can be assumed 
to be finite. Moreover, we can refine {ipk} and assume that supp ipk C U whenever 
supp^ nS 0 ^ 0. Set 

$(z) = Y ^kiz), ip(z) = Y M*) ~ $( z )- 

suppi p k r\S a ^9 k 

The functions ip and ip belong to C'oo(C), satisfy 

suppi/* n S a = 0, suppi/* C U (A. 6 ) 


and 


ip{z) + ip{z) = 1, z £ supp Hi U supp H 2 U supp H 3 - 


(A.7) 


Lemma A.2. Under the conditions of Theorem \ 1.3[ if q is a polynomial, then 
'D^ q (jT) = 0 , where ip is the function constructed above. 


Proof. From (A.7), it follows that q = ipq + ipq on supp/ii U supp/z 2 U supp/ 13 . 


From the definition of T>h in (2.1) we then get 


Using the first condition in (A.6), we get supp(i/*g) D S a = 0, so Lemma 


2.f 


gives 

us 'D^ q (fi) = 0. Since q is a polynomial, we learn from Lemma A.l that T> q (jT) = 0, 
concluding the proof. □ 


We are finally able to prove Theorem |1.3| in its full generality. 


Proof of Theorem \1.3\ Recall the definition of the arcs { 7 ^}, p £ S a , given in (A.4). 
Each of these arcs is a simple contour on the complex plane, and we can find a 
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smooth arc 7 C C for which U 7 p C 7 and C \ 7 is connected. In particular, the 
second condition in (A. 6 ) implies 


supp(^) (~1 (supp (ii U supp H2 U supp ^3) C u Ip C 7. 

p&S a 


(A.8) 


Consider a parametrization 7 : [0,1] —» C of 7 by arc length, set a = 7 ( 0 ) and 
define a continuous function g : 7 — > C by g(z) = 7^7 _1 (z)). That is, </(z) is a unit 
vector tangent to 7 at the point z , varying continuously with 2;. Given ft £ C 2 (C), 
define 

G : 7 —>■ C, G(z) = ^y ( 2 ^) ResO) + ^(z) Im 5(2:)^ 

7 is a simple smooth arc, so it has empty interior. Since C \ 7 is connected 
and G is continuous, Mergelyan’s Theorem tells us that there exists a sequence of 
polynomials ( p n ) converging to G uniformly on 7. In particular, for 7(t) = z this 
implies that the convergence 

Pn 7 Wb (*) = Pn(z)g(z) -t v- Re ff (2:) + — 2:) Img(z) = — (ft(7(t) ) 

oa; 0 ?/ at t= 7 -i( z ) 

holds uniformly on 7, and as a consequence the sequence of polynomials 


5n(A) = / Pn(s)ds + h(a) 

J a 

converges to 

/ -^-(/i(7 (u)))du + ft (a) = h(z) 

J 0 du 

uniformly for z £ 'y. 

In summary, we constructed a sequence of polynomials (q n ) converging uniformly 
to ft on 7 , and for which the sequence of derivatives (q' n ) = (p n ) converges to G 
uniformly on 7 ; in particular there exists M > 0, independent of n, such that 

Wn(x)\ < M, x £ 7 , n> 1 . (A.9) 

Hence, the convergence 

V>( x)q n (x ) -» ip(x)h(x) (A.10) 

holds true uniformly along 7 . Due to ( |X 8 | , this is enough to conclude that the 
convergence above holds uniformly on supp pi U supp H 2 U supp g, 3 , so that 

j ‘b'jiqn d-Hj -£ J dpj, j = 1 , 2 , 3 . (A. 11 ) 


The measures Hi,H 2 and y 3 do not have point masses, so the diagonal {a; = y} 
has zero y,j x ftfc measure. Thus the limit (A. 10 1 also implies that the convergence 


ip(x)q n (x) - ip(y)q n {y) , i>{x)h(x) - ip(x)h(y) 


x-y 


x-y 


holds true pointwise (pij x /afc)-a.e. 

Since the arc 7 is connected, we also know that 


${x)qn(x) - ip{y)qn{y) 


x-y 


< sup 
ae[0,l] 


s((*”)° 7 ) (t) ) 


(A.12) 


t=s 
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whenever i,j/ £ 7 . In virtue of (A.91 and (A.10), the right-hand side in the in¬ 


equality above is uniformly bounded in n. Hence the left-hand side of (A.12) is 


uniformly bounded along 7 , and using (A .81 and once again (A. 10 1 , we can extend 


this conclusion to supp p. 1 U supp ^2 U supp / 13 . Using the Dominated Convergence 


Theorem and (A.12) we conclude 


for j, k = 1,2,3. Combined with ( A.ll|), we finally get 


D $q r (P) V 7,h(P)- 


1ph' 


From Lemma 

ti the other ha 
so from Lemma 


A.2 


we know that (/l) = 0 for every n, hence T>^ h (p,) = 0 . 


ifjh ' 


On the other hand, due to the first condition in (A. 6 ), we have supp(V’Zi) O S a = 0, 
2.8 we get 'D 1 pf l (P) = 0- Thus, 


V h (jl) = 'D$ h {p) + T>^, h (iT) = 0, 


where for the first equality we used (A.7). Since h £ C 2 (C) is arbitrary, Corol¬ 
lary [ 2 ] 2 ] gives us that the measure p is critical, concluding the proof. □ 


Appendix B. Quadratic differentials 

A meromorphic quadratic differential roona Riemann surface 7Z is a differential 
form of type ( 2 , 0 ), given locally by an expression f(z)dz , where / is a meromorphic 
function of a local coordinate z. If z = z(Q is a conformal change of variables, then 

7 (Od( 2 = f(z(0)(dz/d() 2 d( 2 


represents w in the local coordinates (. 

In this Appendix, we sketch the minimal background on quadratic differentials 
used throughout the paper. The general references are the monographs by Strebel 
and Jenkins 
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43 ; some additional information can be found in 67 73 81 


B.l. Critical points and trajectories. The critical or singular points of w = 
fdz 2 are the zeros and poles of /; recall that a zero (resp., a pole) of vo is a point 
p where in a local chart sending p to 0 we have f(z) = z n ip(z ), with ^ 0 , and 
with the integer n > 1 (resp., n < —1). The value n is the order of the critical 
point p, and is denoted by r)(p). The rest of the points of 1Z are called regular, and 
their order is assumed to be rj(p) = 0 . 

Critical points of order < —2 (i.e., poles of order 2 and higher) are called infinite , 
and the rest of the critical points are finite. 

In a neighborhood of any regular point p 1 the primitive 


T (z) 




\/-f(s)ds 


(B.l) 


is well defined by specifying the branch of the square root at p and continuing it 
analytically along the path of integration. Function T (z) provides a distinguished 
or a natural parameter on 1Z in a neighborhood of p. 

We are mostly interested in the trajectories of a quadratic differential w. A 
Jordan arc 7 C 1Z is called an arc of trajectory of tzj if it is locally mapped by T to 
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Figure 22. Structure of trajectories in a neighborhood of a reg¬ 
ular point (left), simple zero (middle) and simple pole (right). 


a vertical line. More precisely, this means that for any point p £ 7 , there exists a 
neighborhood U where the primitive T above is well defined and satisfies 


Re T(z) = const, z € 7 D U. 


(B.2) 


A maximal arc of trajectory is called a trajectory of vj. 

Analogously, the orthogonal trajectories of vj are trajectories of — vj; they can 
be equivalently defined by replacing “Re” by “Im” in (B.2). 


A trajectory 7 extending to a finite critical point along at least one of its di¬ 
rections is called critical ; in the case when it happens in both directions, we call 
this trajectory bounded (also finite or short), and unbounded (or infinite) otherwise. 
Notice that both ends of a short trajectory may coincide, in which case it forms a 
loop on 1Z. 

A w-chain is a connected set on 1Z made of a finite union of arcs of trajectories 
or orthogonal trajectories of w. If no curves in a w-chain belong to orthogonal 
trajectories, we refer to it as a path of trajectories of vj (or a w-path). In this case, 
in order to avoid the trivial situation, two consecutive arcs of a w-path are required 
to intersect at a singular point of vj. 

B.2. The local structure of trajectories. The local behavior of trajectories of 
a meromorphic quadratic differential vj is well understood. 

From a point p of order rj(p) = n > — 1 emanate n + 2 trajectories, forming equal 
angles at p. This covers also regular points, meaning that through any regular 
point passes exactly one trajectory, which is locally an analytic arc (see Figure [ 22 ]). 

An infinite critical point p or order n < — 3 has a neighborhood G with the fol¬ 
lowing property: there are — (n + 2) asymptotic directions, henceforth called critical 
directions , forming equal angles at Pt such that each trajectory entering G 

stays in G and tends to p in one of the critical directions 43 Theorem 3.3]. If 


a trajectory is fully contained in G, then it tends to p in two consecutive critical 
directions (see Figure [23|). 

At a double pole p there are three possibilities. For vj = f(z)dz 2 , we define the 
residue c of vj at z = p to be the residue of \Jf(z) at z = p, which is well defined 
up to a sign. If c € R then there are no trajectories emanating from p and the 
trajectories near p are closed loops. If c £ then there are trajectories emanating 
from c in every direction. In the rest of the cases, the trajectories near p converge 


to p in a spiral form (see Figure 24). 
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Figure 23. Structure of trajectories (solid lines) and critical di¬ 
rections (dashed lines) in a neighborhood of poles of order 3 (left) 
and 5 (right). 



Figure 24. Structure of trajectories in a neighborhood of a double 
pole for c € K (left), c £ zK. (middle) and c £ C \ (R U zR) (right). 


B.3. Global structure of trajectories. There are three possible behaviors for a 
trajectory 7 in the large, 

(i) 7 is a closed curve containing no critical points. 

(ii) 7 is an arc connecting two critical points (which may coincide; in this case 7 
is a closed curve). 

(iii) 7 is an arc that has no limit along at least one of its directions. 


Trajectories satisfying (ii) are called short or finite. Trajectories satisfying (iii) 
are called recurrent , and they are usually a major source of troubles when studying 
the global structure of trajectories of a given quadratic differential. Fortunately 
(for us) in this paper we deal with quadratic differentials with at most 3 poles on a 
genus 0 compact Riemann surface, and the absence of recurrent trajectories in this 
case is assured by Jenkin’s Three Poles Theorem 67 Thm. 8.5, page 226]. 

It is intuitively clear that a non-recurrent and non-closed trajectory 7 has two 
limiting (extremal) values when we travel it in both opposite directions, and that 
can be distinct or equal. For convenience, we will denote these extremal values by 
p(j) and 9 ( 7 ). For instance, for a short trajectory 7 both p( 7 ) and q{ 7 ) are finite 
critical points. 

The set of the critical trajectories of a quadratic differential vj (together with 
their limit points, i.e. the critical points of w) is the critical graph of vj, denoted 
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by Q = Q m . According to 43, Theorem 3.5] (see also 177, §10]), the complement 
of the closure of in TZ consists of a finite number of domains called the domain 
configuration of w. The knowledge of Q (or of the domain configuration of w) is 
sufficient to fully understand the global structure of trajectories of w, as evidenced 
by the following theorem: 


Theorem B.l (Basic Structure Theorem 43 Theorem 3.5]). Let w be a meromor- 
phic quadratic differential on a compact Riemann surface TZ. Suppose in addition 
that w has no recurrent trajectories. 

ThenlZ\Q decomposes into a finite union of disjoint domains UT>, each of them 
bounded by a finite number of critical trajectories. Each domain T> lies into one of 
the following four classes. 

(i) Half plane (or end) domain: It is swept by trajectories converging to a 
pole p of order >3 in its two ends, and along consecutive critical directions. 
Its boundary consists of a w-path with two unbounded critical trajectories and 
a finite number of short trajectories. For some choice of the branch of the 


square root, the natural parameter T in (B.l) is a conformal map from T> to 
a vertical half plane 


H = {w £ C | Rew >c}, 


for some c £ M, and it extends continuously to the boundary of V with the 
identification T (p) = oo. 

(ii) Strip domain; It is swept by trajectories which both ends tend to poles p, q 
of order > 2, possibly with p = q. The boundary dT>\{p, q} is a disjoint union 
of two w-paths, each of them consisting of two unbounded critical trajectories 
converging to p and q, and possibly a finite number of short trajectories. 

For some real constants C\ < c 2 , T maps V conformally to a vertical strip 

S = {w £ C | Ci < Rew < C 2 }, (B.3) 


and it extends continuously to the boundary ofD , with appropriate identifica¬ 
tion of the points p, q with the directions ±ioo. 

(ill) Ring domain; It is swept by closed trajectories. Its boundary consists of two 
connected components, where each of them is a closed w-path. For a suitably 
chosen real constant c and some real numbers 0 < ri < r 2 , the function 
z 1 —> e cT ( z ) maps V conformally to an annulus 

R={w£ C I n < |to| < r 2 }. (B.4) 


and it extends continuously to the boundary ofV. 

(iv) Circle domain; It is swept by closed trajectories and contains exactly one 
double pole, with purely real residue. Its boundary is a closed w-path. 

For a suitably chosen real constant c and some real number r > 0, the 
function z i-A e cT ^ is a conformal map from T> to the circle centered at 
origin and radius r; it extends continuously to dD and sends the double pole 
to the origin w = 0. 


In case w has also recurrent trajectories, a fifth class of domains has to be added 
to the domain configuration of w; we refer the reader to 43 for further details. 
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For a given short trajectory 7 connecting two finite critical points p , q (which 
coincide if 7 is closed) we define its length by 


^( 7 ) = 


r _ 


r _ 

/ \J —vo 

j 7 

= 

Im / \j—vo 

J 7 


> 0 , 


(B.5) 


The width of a strip domain S is defined as a(S) = \c 2 — Ci|, where the constants 
Ci,C 2 are as in (B.3). Alternatively, it can be computed as 


r(S) = 


Re v— vo 


(B. 6 ) 


where we integrate along any path in S connecting two points, p and q , lying on 
different connected components of dS \ {poles of 137}. 

In the same spirit, the width of a ring domain 2? is defined as 


a(V) = 1 log — = 
c r 1 


Re 


/—VO 


where 77 , r^ and c are as in (B.4), and we integrate along any path in 2? connecting 


two points, p and q , lying on different connected components of its boundary. 

For any open simply connected domain 2? C IZ bounded by a tu-chain (that is 
called a vo-polygon ) we define two values: the total order of the singular points of 
vo in 2? is 

Pj ev 

where the summation is along all the singular points Pj of vo in 2?, while the 
contribution from the singular points of w on the boundary is 

«(£>) = P&i)' 

PjGdV 

where the summation is along all the corners p :j of <92?, 

dip) + 2 


Pip) = 1 - OiPj) 


2-7T 


and 9(p) is the inner angle (in radians) at the corner p. 

Both values have a simple relation, as shown by the following simple consequence 


of the argument principle (also known as the Teichmiiller lemma, see 77 Theorem 
14.1]): 


Theorem B.2 (Teichmiiller lemma). If V is a vo-polygon, then 

k{ V) = 2 + \{V). 


(B.7) 


A straightforward corollary of this lemma (and also a direct consequence of the 
maximum principle for holomorphic functions) is the following fact: 

Corollary B.3. If vo is analytic (has no poles) in a vo-polygon T>, then <92? must 
contain at least one pole of vo. 


This corollary is also the basis for the general principle P.3 in Section [4. 5. 1| 
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Figure 25. Numerical evaluation of the critical graph for r = 1/25. 


Appendix C. Numerical experiments 


In this Appendix we look under the hood of the calculation of the functions uij in 


(4.19), as well as of the numerical procedures used to find and plot the trajectories 
of the quadratic differential w from Section [4] 

As it was mentioned, Figure[4]was obtained by means of a numerical evaluation of 
the integrals defining the functions ujj in ( |4.19[ ). For that, we compute the integrals 
of the form 

r z 2 

/ Zj(s)ds 


along the line segment joining chosen roots Z\. z 2 of the discriminant in (3.5) by 
means of the composite trapezoidal rule. 

For r fixed, we compute the points Zk = Zk{r), choose a value m £ N, and 
consider a grid of m + 1 equally spaced nodes {pj I”!/} 1 , 


j — 1 

Pj = ~i H- (z 2 -zi), i = !,••• ,to+1, 


At 


that will be used as the quadrature nodes for the composite trapezoidal rule, 
each z = pj we solve (3.9) numerically, obtaining an (unordered) set of three 
solutions Cfe(Pj)- Comparing their values at consecutive quadrature points, we 
collect them into a sequence of (ordered) vectors Vj = (Pj)i 

j = 1,..., m + 1, where a is a permutation of {1,2,3} that does not depend on j. 
We can then determine the permutation a using equations (4.8) and (4.14) (4.17) 
as boundary conditions. 

This procedure is repeated for n equally spaced values of r, 

k 


T = T k = 


4(n +1) ’ 


k = 1, 


(notice that £ [0,1/4); we stop at k = n in order to avoid the degenerate situation 
at r = 1/4 for which some of the endpoints of integration diverge to oo). 

The result of these calculations with n = 1000 and m = 10000 is plotted in 
Figure [4] 

On the other hand, we also performed numerical experiments that helped us to 
build the intuition to predict (and confirm) the structure of the critical graphs in 
Section [4j A sample of such graphs is presented in Figures [25] [27] 
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Figure 26. Numerical evaluation of the critical graph for r = 1/5. 



Hi 




Hi 


Figure 27. Numerical evaluation of the critical graph for t = 1/5. 


The numerical procedure used for these pictures is as follows. If 7 (t),t £ JcK, 
is a parameterization of a trajectory of w = Q 2 (z)dz 2 , then 

Q 2 ('v(t)W(t) = if(t), teJ, 

where / is some real valued function. A different choice of / correspond just to a 
reparameterization of 7 . In a natural (arc-length) parametrization of 7 the equation 
above takes the form 


7 '(i) = 


Q 2 (7W) 


f G J. 


(C.l) 


IQ 2 (7(*))r 

This is an ordinary differential equation for 7 that can be solved by standard 
methods (e.g. the family of Runge-Kutta algorithms). 

We should bewared of two aspects when implementing this method. The first one 
concerns the initial values, that usually are at branch points of (3.91, and thus, at 


zeros of Q , from where more than one trajectory emanates. This can be resolved by 
perturbing the initial point along the prescribed direction, which can be obtained 
from a local analysis of Q at the singular point. 

The second aspect concerns the choice of branches of Q 2 in (C.l), where we 


take advantage of the fact that all the branch points of (3.9) are quadratic. The 
value of Q 2 is found numerically by solving (13.91), which gives us the three possible 


values Near a branch point two of these values are close (corresponding to the 
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two solutions coinciding at this branch point), allowing us to distinguish the third 
solution. This, in turn, allows to recognize on the the possible branches of Q 2 , but 
not the other two. 

For instance, at 2 = &2 we will have £1 ~ £3 (see pT 7 | l), which singles out the 
branch £2 and makes the branch Q 2 — (£1 — £ 3) 2 easily distinguishable. In order 
to identify the remaining two branches (i.e., (£2 — £i ) 2 and (£2 — £a) 2 ) we must 
use further results about the local and global structure of the critical graph of w, 
established in Section [4j 
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